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ABSTRACT 


Experimental and theoretical research is being undertaken on 2-D 
ax i symmetric geometries under low speed, nonreacting, turbulent, swirling flow 
conditions. The flow enters the test section and proceeds into a larger 
chamber (the expansion ratio D/d = 2) via a sudden or gradual expansion (side- 
wall angle a = 90 and 45 degrees). Inlet swirl vanes are adjustable to a 
variety of vane angles with values of <t> = 0, 38, 45, 60 and 70 degrees being 
emphasized. 

This research is being funded by NASA Lewis Research Center and AF Wright 
Aeronautical Laboratories under Grant Number NAG 3-74 since July 1, 1980. The 
present document is the Fourth Semi-Annual Status Report and discussion is on 
those activities undertaken during the second half of the second year of the 
study. Section 1 of this document describes the recent activity, while 
Section 2 gives a brief summary of the present status and a list of recent 
publications. Appendices A through E contain an M.S. Thesis, two recent 
internal reports, and two recent AIAA papers. 
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1. RECENT PROGRESS 


New features of the present year's study are: the turbulence 

measurements are being performed on swirling as well as nonswirling flow; and 
all measurements and computations are also being performed on a confined jet 
flowfield with realistic downstream blockage. Recent activity in the research 
program falls into three categories: 

1. Time-mean flowfield characterization by five-hole pitot probe 
measurements and by flow visualization. 

2. Turbulence measurements by a variety of single- and multi -wire hot- 
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wire probe techniques. 

3. Flowfield computations using the computer code developed during the 
previous year's research program. 

During the last six months, students have drawn up six brief papers 
covering our studies. They were supplied to the Contract Monitor with the 
April 1, 1982 newsletter. Three internal documents relating to degree 
requirements were successfully written, presented and approved by thesis 
Committees.^" 9 Hyung K. Yoon recently completed his M.S. Thesis, included 
here as Appendix A. Predictions resulting from very accurate inlet flow 
measurements 5 (see Appendix B) are described and discussed in Ref. 11 (see 
Appendix C). Four ASME or AIAA conference papers were prepared for June 1982 

meetings, the two most relevant of which are included in Appendices D and 
E.12,13 


1.1 Time-Mean Velocity Measurements 


A recently completed M.S. Thesis* 0 investigated nonswirling and swirling 
flows (reaching higher swirl strengths than previously) in the test section. 
Later measurements included the effect of contraction nozzles of area ratio 2 
and 4, located at two axial positions. The effects of several geometric 


parameters on the flowfield were investigated including: side-wall expansion 
angle a - 90 and 45 degrees, swirl vane angle 4* = 0, 38, 45, 60, and 70 
degrees, and contraction nozzle location L/D = 1 and 2 (if present). Data 
acquisition was via a five-hole pitot probe enabling three time-mean velocity 
components in the axial, radial and azimuthal directions to be measured. 
Velocities, yaw and pitch angles are extensively tabulated and plotted in a 
directly usable form. Artistic impressions of recirculation zones are then 
deducible. Full details are included in the M.S. Thesis included here as 
Appendix A. 

Some results are now described briefly. The nonswirling confined jet 
possesses a corner recirculation zone extending to just beyond x/D = 2 with no 
central recirculation zone. The presence of a swirler shortens the corner 
recirculation zone and generates a central recirculation zone followed by a 
precessing VGrtex core. The effect of a gradual inlet expansion is to 
encourage the flow to remain close to the side wall and shorten the extent of 
the corner recirculation zone in all cases investigated. A contraction nozzle 
of area ratio 2 has little effect on weakly swirling and strongly swirling 
flows, which are dominated by forward flow and centrifugal forces, 
respectively. For the intermediate swirl case of = 45 degrees, it 
encourages forward movement of otherwise slow-moving air and thereby shortens 
the central recirculating zone. A strong contraction nozzle of area ratio 4 
has a more dramatic effect on the flowfields. 

1.2 Swirler Performance 

Five-hole pitot probe measurements have also been taken immediately 
downstream of the swirler, with expansion blocks removed, in contrast to 
earlier measurements which were made at the expansion plane located 3.2 cm 
further downstream. This study defines better the performance characteristics 
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of the swirler. Measurements of time-mean axial, radial and tangential 
velocities at the swirler exit plane for a range of blade swirl angles <J> from 
0 to 70 degrees show clearly the effects of centrifugal forces, recirculation 
zones, and blade stall on the combustor inlet velocity profiles. A brief 
paper^ describing these activities is included in Appendix B; more complete 
documentation is currently being drawn up. Results show that exit velocities 
are functions of radial and azimuthal parameters. 

Assumptions of flat axial and swirl profiles with radial velocity equal 
to zero are progressively less realistic as the swirler blade angle 
increases. At moderate swirl <f> = 45 degrees, linearly increasing profiles of 
u and w with radius are appropriate with v nonzero. At strong swirl <P = 70 
degrees, even more spiked profiles are appropriate with most of the flow 
leaving the swirler near its outer edge. Non ax i symmetry prevails in all swirl 
cases investigated. 

1.3 Computational Studies 

STARPIC is an advanced computer code developed earlier at Oklahoma State 
University to predict corresponding confined swirling flows to those studied 
experimentally. The calculation method includes a stairstep boundary 
representation of the expansion flow and a conventional k- e turbulence 
model. The predictions include velocity profiles, recirculation zone 

characterization and mean streamline patterns. Currently, predictions are 
being made with the standard k-e model and no downstream blockage. 
Predictions now include realistic inlet conditions for a complete range of 
swirl strengths (as described in Appendix B and other unpublished data), and 
results conform more realistically to measured flowfield velocities and 
streamline patterns. A full report^ on these activities appears as Appendix C. 
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The appended report provides numerical predictions using various inlet 
velocity starting conditions for the case of swirl vane angle equal to 45 
degrees. The validity of flowfield predictions resulting from the choice of 
inlet profiles is assessed by comparing the predicted velocity profiles with 
corresponding experimental velocity profiles. Results demonstrate that 
realistic predictions are forthcoming only from the inclusion of realistic 
axial, radial and swirl velocity profiles as inlet conditions, and that 
considerable errors occur if unrealistic idealized inlet conditions are 
used. Predictions are then exhibited for a range of swirl strengths (<j> = 0, 
38, 45, 60 and 70 degrees) using measured inlet axial, radial and swirl 

4 

velocity profiles in each case. The ensuing flowfields are characterized via 
velocity profiles and streamline patterns, and illustrate the large-scale 
effects of inlet swirl on flowfields. 

Predictions of downstream blockage effects (two blockage sizes at two 
axial locations) have now been completed and these are being assessed in the 
light of the experimental data given in Appendix A. The amalgamation of a 
full report on these activities is currently underway. 

Oblique to the main thrust of the present research program is application 
of the STARPIC code to MHD combustor flows, simple numerical solutions to the 
bluff-body flameholder wake problem, and prediction of 3-D turbulent flows. A 
brief research paper ^ on the second of these is included as Appendix D. 
There, numerical predictions are made of recirculation zones behind bluff -body 
flame stabilizers, showing quantitatively the effects of forebody geometry, 
blockage ratio, lateral position of the blockage and inlet swirl on the 
central recirculation zone. A simple transient Navier-Stokes solution 
algorithms and laminar flow simulation are used with 'free slip' and 'no slip' 
wall boundary conditions, thus illustrating how a basic approach may be used 
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to solve a sophisticated fluid dynamic problem. A 3-D turbulent version of 
this code is described in the research paper of Appendix D of the Third Semi- 
Annual Status Report. Turbulence is simulated by means the two-equation k-c 
turbulence model; species diffusion and buoyancy are also included. Two 
applications of the code are presented to local destratification near the 
release structure of a reservoir aid to the deflection of a Jet entering 
normally into a uniform cross-flow. Predicted results exhibit good agreement 
with experimental data, showing that a useful characterization of fully three- 
dimensional flows is now available. 

1 . 4 Hot-Mire Anemometry 

Major attention has been given to application of the one-wire hot-wire 
techniquein the present confined jet facility. Determination of all time-mean 
vleocity components and all normal and shear Reynolds streses results from the 
6-orientation mean and fluctuating data at each f lowfield location. The 
experiments are being performed to provide the information necessary for 
turbulence modeling development in the confined jet facility. Recently 

swirling flows have been investigated (with <f> = 38 and 45 degrees) and a 
research paper* 3 on the results, including the earlier nonswirling results, is 
included as Appendix E. Excellent cost-effective results are presented, with 
comparisons with independent data illustrating the reliability of the 
technique. Finally, a sensitivity analysis of the data reduction technique 
has been undertaken; it forms the major ingredient of an uncertainty analysis. 

A new cross-wire probe is currently being used to make a more accurate 
evaluation of the turbulent shear stress than previously possible. This 
technique is being used for the nonswirling flow without downstream blockage, 
and is soon to be used with downstream blockage. These data will enable 
turbulent viscosity components to be deduced (from the turbulent stress vs. 
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mean velocity spatial gradients) and will aid the turbulence modeling aspects 
of the study. 

A series of hot-wire measurements have been planned, using a three-wire, 
hot-wire probe with direct computer interface and data reduction. These 
measurements are to be performed by subcontractor Dynamics Technology, Inc., 
under the direction of Dr. Dennis K. McLaughlin. There, experiments are 
progressing behind schedule and recently a no-cost time extension to October 
22, 1982 was agreed between them, Oklahoma State University and NASA Lewis 
Research Center. For this reason, no report on their activity is included in 
the present document, but copies will be distributed as. soon as they are 
available. 

1.5 Construction Activity 

Dynamics Technology has been assembling for their subcontract and other 
activity a swirling jet facility which will duplicate the flowfield produced 
by the Oklahoma State University facility. However, it will be instrumented 
differently, with the three-wire hot-wire probe and its three-degree of 
freedom probe drive being major constituents. The complex probe drive, swirl 
pack and expansion blocks have been built by students at Oklahoma State 
University and loaned to the subcontractor for their experimentation. 
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2 . 


PRESENT STATUS 


2.1 Summary 

During the second year of the present study, new features include: 

X. Turbulence measurements are being performed on swirling as well as 
nonswirling flows. 

2. The effect of size and location of a downstream blockage is being 
investigated. 

Successful studies include: 

1. Flow visualization. 

2. Five-hole pitot probe measurements. 

3. Single-wire and multi-wire hot-wire measurements. 

4. Flowfield computations from realistic inlet conditions. 

Proposed future studies include: 

X. Expansion Ratio D/d = X.5 - in addition to the D/d = 2 value 

currently being used, including flow visualization, time-mean 
velocity and normal and shear turbulent Reynolds stress measurements. 

2. Corresponding Predictions - using advanced turbulence modeling 

developments. 

3. Inlet Turbulence Level - effects to be clarified via flow 

visualization, and some time-mean and turbulent stress measurements, 
for a swirling flowfield. 

4. Turbulence Length Scale - to be measured with one and/or two hot-wire 
probes, in nonswirling and one swirling flow. 

5. Turbulence Model Developments - in the light of on-going activity 

both at Oklahoma State University and elsewhere. 

The completion of these tasks, and the accomplishment of the goals and 

objectives of the grant will make a substantial contribution to the 

understanding and prediction capability of complex turbulent flows. 
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2.2 Publications 


The following list of publications covers research conducted with the 
support of the grant, and does not include those listed in previous Semi- 

Annual Reports: 

1. Abujelala, M. T., Confined Swirling Flow Predictions: Effect of 

Inlet Velocity Assumption. Proceedings of 13th Southwestern Graduate 
Research Conference in Applied Mechanics, held at the University of 
Oklahoma, Norman, OK, on April 16-17, 1982, pp. 267-273. 

2. Busnaina, A. A., Turbulent Three-Dimensional Flowfield Predictions. 

Ibid, pp. 17-23. 

3. Jackson, T. W., A Single Hot-Wire Technique for Swirling Flows, 

Ibid, pp. 260-266, 

4. McKillop, B. E., Crossed Hot-Wire Anemometry in Turbulent Flows. 

Ibid, pp. 11-16. 

5. Sander, G. F., Annular Vane Swirler Performance. Ibid, pp. 274-280. 

6. Yoon, H. K., Five-Hole Pitot Probe Time-Mean Velocity Measurements in 
a Combustor Flowfield. Ibid, pp. 131-137. 

7. Abujelala, M. T.» Turbulence Models and Swirling Flow Prediction. 

Ph.D. Preliminary Write-Up, MAE, Oklahoma State University, June 
1982. 

8. Jackson, T. W., Turbulence Measurements in Axisymmetric Combustor 

Geometries. Ph.D. Preliminary Write-Up, MAE, Oklahoma State 
University, June 1982. 

9. Busnaina, A. A., Transient Three-Dimensional Predictions of Turbulent 

Flows in Cylindrical and Cartesian Coordinate Systems. Ph.D. 

Qualifying Write-Up, MAE, Oklahoma State University, June 1982. 
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10. Yoon, H. K., Five-Hole Pitot-Probe Time-Mean Velocity Measurements in 
Confined Swirling Flows, M.S. Thesis, MAE, Oklahoma State University, 
July 1982. 

11. Abujelala, M. T, and Li 1 ley, D. 6., Confined Swirling Flow 

Predictions: Effect of Swirl and Inlet Profile Assumptions, Report, 

MAE, Oklahoma State University, June 1982. 

12. Vatistas, G, H., Lin, S., Kwok, C. K. and Lilley, D. G., Bluff-Body 

Flameholder Wakes: A Simple Numerical Solution. Paper AIAA-82-1177, 

Cleveland, OH, June 21-23, 1982. 

13. Janjua, S. I., McLaughlin, D. K., Jackson, T. W. and Lilley, D. G., 

«r 

Turbulence Measurements in a Confined Jet Using a Six-Orientation 
Hot-Wire Probe Technique. Paper AIAA 82-1262, Cleveland, OH, June 
21-23, 1982. 
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CHAPTER I 


INTRODUCTION 

1.1 The Investigation of the Flowfield in 
a Confined Jet 

A need for more thorough understanding of the fluid dynamics of the 
flow in a gas turbine combustor chamber has been recognized by designers 
in recent years. The flowfield is presently being investigated using 
various methods of approach, such as computer simulation (1-5) , .flow- 
visualization (5-6) and time-mean velocity measurements (7-14) for both 
swirling and nonswirling flows. 

Until recently, extensive and reliable experimental data for swirl- 
ing flow have been unavailable, and more experimental data are needed. 

Such data together with combustion characteristics would be very useful 
to designers as well as to theoretical investigators who could compare 
the experimental data with the data obtained from the simulated mathemati- 
cal methods. 

Current research focuses on measuring time-mean velocity and investi- 
gating the extent of the recirculation zone (known to have great effects 
on flame stability and length) and the effects of geometric parameters on 
the flowfield in a confined can-type gas turbine combustor. 

In order to study the flow phenomena in the combustor, it is neces- 
sary not only to know the static pressure and flow velocity, but also the 
flow direction. In addition, the sensing probe must be sufficiently 
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small to be capable of point measurement without appreciably affecting 
the flow in the passages by its presence. In these experiments, the five- 
hole pitot probe was used because it is acutely sensitive to flow, con- 
venient to manipulate and relatively compact. Measurements were taken 
for nonsirling and swirling flow with and without a contraction block in 
an axlsymmetrlcal confined jet facility simulating typical geometry. All 
flowfields were nonreacting. 

1.2 Review of Previous Flowfield 
Measurement Studies 

Lee and Ash (7) developed and calibrated a three-dimensional spher- 
ical pitot probe which would measure the static pressure and the magni- 
tude and direction of the velocity vector for any arbitrary flow angle 
without needing to be adjusted. It was designed mainly for measuring the 
relative flow through rotating three-dimensional blade cascades. 

An investigation of flow patterns in the gas turbine combustion 
chamber with a three-dimensional “robe was done by Hiett and Powell (8) 
in the early nineteen sixties. They found that there is a large recir- 
culation zone between the main flow along the centerline and the second 
flow along the chamber wall, and that there is a high temperature zone 
in this area. 

Janjua (9) measured three velocity components and turbulence levels 
for nonswirling flows with the six-orientation hot-wire technique. An 
uncertainty analysis done on this technique reveals that turbulence in- 
tensities and shear stresses are extremely sensitive to some input par- 
ameters such as yaw factor and mean voltages. 

Velocity distribution in a sudden expanding circular duct was 
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investigated by Moon and Rudinger (10) with a Laser-Doppler velocimeter. 
These experiments located both the dividing streamline, which separates 
the main flow from the recirculation region, and the point of the flow 
reversal In the recirculation region. 

Chaturvedi (11) measured the mean velocity and pressure and the char- 
acteristics of the turbulence by using a stagnation tube and pitot probe 
in a can-type axisymmetric expansion combustor for isothermal nonswirling 
flowfields. His data will be compared with the current research results 
because of the close geometric similarity and experimental conditions. 
Furthermore, in his experiment a single and cross-hot wire were used to 
measure the longitudinal component of turbulence intensity, u * . radial 
intensity, y* , and turbulent shear stress, u’v* . 

Systematic research of swirl was done by Mathur and MacCallum (12). 
The test facilities which they used for the free jet consisted of two 
concentric tubes with supplies of air to both the central (primary) and 
annular (secondary) streams. They used a three-dimensional probe to 
measure the velocity components in the three-dimensional flow. They 
measured the pressure drop across the vane swirler by both experimental 
work and theoretical analysis. They also investigated the distribution 
of the static pressure along the axis of the jets and the variations of 
the velocity distributions in the tube according to the change of swirl 
vane angle in the free jet. Mathur and MacCallum (12) followed the same 
procedure in the experiments for the enclosed jets as for the free jets 
and observed similar physical characteristics for the flows according to 
geometric changes. They determined that a large central recirculation 
exists along the axis in the swirling free and confined jets. A combus- 
tion chamber of square cross-section was used. 
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Pratte and Keffer (14) have Investigated a swirling jet having a 
moderate ratio of swirling to axial momentum. They showed that the flow 
achieved a self similarity for the mean velocities rather quickly while 
the normal turbulent intensities reached a self-similar state after a 
longer period of jet development. Also the entrainment rate and angle 
of spread for the swirling jet was found to be nearly twice that of the 
nonswirling free jet. 

1,3 Flowfield Independence of Reynolds Number 

Although the current experiments are being done under low speed and 
nonreacting conditions because of the limitations of the experimental 
facilities, the results can be applied to actual combustor hardware, 
which usually operate at higher Reynolds numbers than those used in the 
present test facility. 

Moon and Rudlnger (10) noticed that the location of the flow rever- 
sal is practically independent of the flow velocity over the range from 
about 40 to 90 m/s in the small tube. They also observed that, as the 
Reynolds number increases, the uncertainty in determining the reattach- 
ment point appears to decrease. 

Rhode (5 ) , who performed experiments to investigate the asymptotic 
invariance of Re^ on the flowfield, showed that there is little change 
in the flowfield for Re^ > 47,000, where Re^ is based on the nozzle inlet 
velocity and diameter. His study was done for x/D * 2.5, a * 90 degrees, 
and iji » 38 degrees. Because the flow was considered to be most sensi- 
tive to Re d at a location slightly downstream of the recirculation zone 
and the recirculation zone was observed just before the axial station of 
2.5. Also the central rather than corner recirculation zone was selected 
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since most of the combustion and flame stability occur there (17-19) » 
Measurement traverses wore repeated for the same conditions while increas- 
ing Re^ for each of the f lowf ields . in the present study the experiments 
were repeated using the same procedure as Rhode but for x/D * 0.5, ct * 

90 degrees, and <}» * 45 degrees, where the most extensive and strongest 
turbulent reverse flows had been observed in previous experiments. The 
results indicate also the excellent consistency of velocity profiles in- 
dependent of Reynolds numbers for Re^ > 90,000 at the axial location as 
shown later in Figure 11. The present study uses Reynolds numbers rang- 
ing from 53,000 to 220,000 as tabulated in Table X. For strong swirling 
fiow with swirl vane angle <}> * 70 degrees, the Reynolds number is 53,000. 
This is not high enough to be sure that the flowfiled is independent of 
Reynolds number, but 53,000 is the highest value which can be reached 
for t|) « 70 degrees in the current experimental facilities. In this study, 
a flowfield is assumed to be Independent of Reynolds number for Re^ > 
53,000. 


1.4 The Scope and Significance of the 
Present Research 

In the present study, a five-hole pitot probe was used to measure 
the magnitude and direction of the time-mean velocity in the combustor. 
Measurements have been carried out for sudden expansion and gradual ex- 
pansion (side-wall expansion angle a ** 90 and 45 degrees) for nonswirl- 
ing and swirling flows [swirl vane angle zero (swlrlet removed), 38, 45, 
60, and 70 degrees];; and then with contraction block inserted at certain 
axial stations (L/D * 1.0 and 2.0). The raw data measured by a five-hole 
pitot probe are reduced by a FORTRAN computer program (which is described 
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in Appendix G) and plotted as axial, radial, and swirl velocity compon- 
ents by a plotting program using a complot plotter. 


The current research is significant for the following reasons (23) : 

1. The magnitude and direction of the flow velocity is measured. 

2. The extent of the corner and central recirculation zones is 
identified. 

3. The effects of geometric parameters on the flowfields in the 
combustor chamber are determined, 

4. Miscellaneous results such as the performance of the swirl vanes 
and the independence of the flowfields from the Reynolds number are pro- 
vided. 

The research is limited by the following factors: 

1. The difficulty in investigating the exact reattachment points 
of the flows. 

2. The low sensitivity of the probe to small velocities. 

1.5 Outline of the Thesis 

In the previous sections, the significance and scope of this study 
are introduced. Basic concepts are also given with reference to previous 
studies . 

Chapter II describes the facilities and instrumentation employed. 
This chapter also explains the operation of the five-hole pitot probe, 
thCi concept of swirling flow, and the way a swirler performs. 

The procedures of data measurements and reduction are described in 
Chapter III. A separate documentation in Appendix C contains a detailed 
description of the data reduction computer program and instructions for 
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its use. Descriptions of the calibrations of the voltmeter and the five- 
hole pitot probe are also included in this chapter. 

Results for sudden expansion flows are discussed with emphasis on 
the extent of the recirculation zone in Section 4,1. Effects of geomet- 
ric parameters such as a gradual expansion and a contraction nozzle are 
described in Sections 4.2 and 4.3. 

Chapter V summarizes all of these results and suggests further re- 
search, 

Appendix A and B include tables and figures, respectively. User’s 
guide to five-hole pitot probe data reduction computer program is in 
Appendix C. 



CHAPTER II 


EXPERIMENTAL FACILITIES AND FIVE-HOLE 
PITOT PROBE INSTRUMENTATION 

2.1 An Idealized Flowfield 

The test facility used to obtain a uniform flow of relatively low 
turbulence level is a wind tunnel designed and built at Oklahoma State 
University. The test section of the swirling confined-jet facility is 
a simulation of a typical axisymmetric can-type combustion chamber from 
a gas turbine engine as shown in Figure 1. A schemacxo of the facility 
is given in Figure 2. Since the actual gas turbine combustor is too com- 
plicated to be simulated in a laboratory, it is simplified, but the typ- 
ical geometric features such as the presence of a sudden or gradual 
expandion block, a contraction nozzle, a swirler and a can-type shape of 
the flowfield domain are kept, 

Ambient air enters the low-speed wind tunnel through an air filter 
consisting of a large box covered with foam rubber. Next a six-blade 
propeller fan, driven by a 5 h.p. varldrive motor which can be varied con- 
tinuously from 1600 rpm to 3100 rpm, accelerates the air to the desired 
velocity. 

The intensity of the turbulence caused by the driving fan is partial- 
ly reduced by gradually expanding the flow into a stilling chamber con- 
taining numerous fine mesh screens. The turbulence intensity is further 
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reduced by passage through a turbulence management section. This section 
consists of a perforated aluminium plate (2 mm diameter holes) followed 
by a fine mesh screen, a 12.7 cm length of packed straws, and five more 
fine mesh screens. Most of the control of the turbulence occurs in this 
portion. The airflow, which is now relatively uniform and of low turbu- 
lence, then enters a contoured nozzle, which was designed to produce a 
minimum adverse pressure gradient on the boundary layer and thus avoid 
local separation and flow unsteadiness. The ratio of the area of the 
maximum cross section of the contour nozzle to that of the nozzle throat 
is approximately 22.5. The diameter of the nozzle throat is about 15 cm. 
There is a hole of 1 cm diameter just upstream of the nozzle throat for 
insertion of a standard pitot-static probe, so that the dynamic pressure 
can be measured upstream of the swirler. With the measured dynamic pres- 
sure, the corresponding nozzle throat velocity and Reynolds number are 
obtained from the dynamic pressure conversion characteristic shown in 
Figure 3. 

2.2 Test Section 

The test section is composed of a swirl vane assembly, an expansion 
or contraction block, a contraction nozzle, and a simplified can-type gas 
turbine combustion chamber. Details concerning the swirl vane assembly 
are given in Section 2.4. 

The expansion or contraction block can be used optionally to give 
the simulated geometric characteristics desired for a certain expansion 
condition. The expansion blocks are interchangeable and those currently 
in use have expansion angles of 90 and 45 degrees. The expansion blocks 
have an inside diameter the same as that of the nozzle throat, and an 
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outside diameter approximately twice that of the nozzle throat, to match 
the test section inside diameter. The experimental results of Rhode (5) 
show that varying the expansion angle of the side wall does not affect 
the flowfield much except that the corner recirculation region decreases 
sharply as the expansion angle is decreased. The current study was con- 
ducted with side-wall expansion angles of a * 90 and 45 degrees, to in- 
vestigate the effects of both sudden and gradual expansions. A contrac- 
tion block is used to simulate the confining wall of a real gas turbine 
combustor at the downstream end. The block was designed s#ith a rounded 
profile that is one-quarter of a circle and a downstream cross-sectional 
area half that of the test section. The rounded upstream face is smooth- 
ly finished to reduce the surface shear stress. Details of the contrac- 
tion block employed are given in Figure 4. 

The test section used in this study is a straight, long circular 
plexiglass tube having a 30 cm inside diameter. Holes with 2.54 cm di- 
ameter are located along the length of the tube at distances which are 
multiples of x/D - 0.5, so that the pitot probe can be inserted into the 
flow and traverses made in the radial direction. It does not have film 
cooling holes or dilution air holes. The test section is carefully 
aligned with a low-power laser beam so that the centerlines of the test 
section and the wind tunnel coincide. 

2.3 Five-Hole Pitot Probe Instrumentation 

A five-hole pitot probe is one of the few instruments capable of 
measuring both the magnitude and the direction of fluid velocity simul- 
taneously. The five-hole pitot probe used In this study is a model 
DC-125-12-CD manufactured by United Sensor and Control Corp. The 
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sensing head is hook-shaped to allow for probe shaft rotation without 
altering the probe tip location. Little information is available con- 
cerning the effects of turbulence on a pressure probe in a swirling flow. 
However, it is asserted that the five-hole pitot probe is accurate with- 
in approximately 5 percent for most of the measurements (5) . This value 
may increase to 10 percent as the velocity magnitude falls below approxi- 
mately 2.0 m/s because of the Insensitivity of the probe to low dynamic 
pressure and the dependence of probe calibration on the probe Reynolds 
number. 

There are three standard methods of operating the five-hole pitot 
probe: 

1. To adjust the orientation of the probe in both pitch and yaw so 
that the probe is aligned with the local flow direction. 

2. To determine the flow direction from the calibration relation- 
ship between probe pressures and flow direction while maintaining a fixed 
probe orientation. 

3. To align the probe yaw angle with the flow yaw angle while de- 
ducing the pitch angle and total velocity coefficient calibration 
characteristics . 

The third method was employed in this study because it used readily 
available orientation equipment and relatively simple data reduction 
procedures . 

The instrumentation assembly, in addition to the five-hole pitot 
probe, is composed of a manual traverse mechanism, two five-way ball 
valves, a differential pressure transducer, a power supply, and an inte- 
grating digital voltmeter. The differential pressure transducer is model 
590D from Datametrics, Inc. It has a differential pressure range of from 
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0 to 1.3 x 10 N/to , The integrating voltmeter is the TSI model 1076. 

As auxiliary equipment, a model 631-B strobotac from General Radio, Inc. 
was used to check the fan speed and a micro-manometer along with a pitot- 
static probe was used to measure the dynamic pressure in the nozzle 
throat just upstream of the swirler . Also, a barometer/thermometer unit 
from Cenco Corporation was used for local pressure and temperature read- 
ings. Measurement apparatus are shown in Figure 5. 

2.4 Swirler 


Swirling jets are used as a means of controlling the length and 
stability of the flame in many combustor designs. The effects of the 
degree of swirl on the flowfields and combustion have been extensively 
investigated (5, 12, 19-21). A swirl number characterizes the degree of 
swirl. It is usually defined as the ratio of axial flux of tangential 
momentum to the axial flux of axial momentum (including a pressure con- 
tribution) , nondimens ionalized by use of the inlet nozzle radius. Thus; 


S 


C x d/2 


(2-1) 


where: 

G 0 * the axial flux of tangential momentum flux; and 
Sc * t ^ ie ax -*- a ^ fl ux of axial momentum, including pressure contri- 
bution. 

Often the pressure term is neglected, and the expression reduces to (20) 


„ 2 ,1 - Z~\ . 

S * -z ( r) tan <{> 

J 1 - Z l 


( 2 . 2 ) 


where Z 


d h /d for the case of an annular vane swirler idealized with 
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flat axial and flat swirl velocity profiles at the inlet. This is the 
case in the present study. 

The swirler used currently was constructed at Oklahoma State Univer- 
sity according to preliminary design plans from NASA, The schematic of 
the swirler is shown in Figure 6. It has ten adjustable vanes and a hub 
with a streamlined nose and a flat downstream face. The nose has a 
hyperbolic shape with a very smooth surface to reduce the pressure drop. 
The flat face and location of the hub is a simulation of an actual fuel 
injector. The flat blades are wedge-shaped to give a constant pitch-to- 
chord ratio of approximately one which should give good turning effici- 
ency (21) . 

Mathur and MaeCallum (12) thought that an overlap between the adja- 
cent vanes is a determinant parameter to give a complete deflection of 
air. They found that a positive overlap is essential for effective 
directing of the air and a 30 degrees overlap is adequate for complete 
deflection. An approximately variable 20 degrees overlap is provided by 
the swirler used in this study. The positive axial velocity which was 
observed near to the hub might be the result of an insufficient overlap. 
An extensive study concerning the swirler performance is being conducted 
by Sander (24) . 
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CHAPTER III 

MEASUREMENT AND REDUCTION OP DATA 

A schematic of the five-hole pitot probe is shown in Figure 7 , To 
obtain all the components of the velocities in the axial, radial, and 
azimuthal directions at a given location in the model combustor, one must 
know the yaw angle g, the pitch angle 6, and the magnitude of the total 
velocity V at that point. The yaw angle (3 can be measured directly from 
the reading of the rotary vernier. The pitch angle 6 and the total velo- 
city V can be found from the corresponding calibration in a free jet. 

With the above information, Rhode (5) wrote a data reduction comput- 
er program. This program reduces the raw pressure data from direct 
measurement into u, v, and w velocity component data. 

3.1 Measurement Procedures 

When using the five-hole pitot probe, one aligns the probe tip with 
the direction of the air flow in the horizontal plane and then reads the 
yaw angle g and the two differential pressures (P N - P g ) and (P c - P w ) . 
These pressures are used to obtain the pitch angle 6 in the vertical 
plane and the magnitude of total velocity vector from the calibration 
characteristics . 

The first step in preparing to take measurements with the five-hole 
pitot probe is to measure the dynamic pressure of flow with a pitot- 
static probe at the nozzle throat just **cforo the sv/lrler. Then the 



14 


15 


pressure transducer must be nulled for zero differential pressure. This 
is done by turning an adjustment screw on the transducer until the volt- 
meter reading for zero differential pressure (both ports open to atmo- 
sphere) is zero. 

The next step in preparing for measurements is to adjust the posi- 
tion of the test section until the wind tunnel and test section axes 
coincide, and also to adjust the traverse mount until the probe zero 
point lies on the test section axis. Then the varidrive motor is set at 
the proper rpm to give the desired air velocity at the nozzle throat. 

The fan speed is checked regularly with a strobotac to keep the air ve- 
locity constant. 

Before the production measurements are taken, the five-hole pitot 
probe rotary vernier must be set for zero yaw angle so that the x and 
0 axes of the measurement coordinate frame coincide with those of the 
test section. This is achieved by adjusting the yaw angle to be zero at 
the center of the test section inlet for nonswirling flow. Each radial 
traverse measurement begins at the centerline, which is identified with 
a low-power laser beam. The pressure and yaw angle data are read at 0.3 
inch increments up to 5.7 inches from the centerline. In this study, 10 
seconds time-mean data readings are used. To measure 10 seconds time- 
mean values, 30 to 40 seconds settling time is needed to allow transients 
to subside in the measuring apparatus. 

The first measurement for each location is the yaw angle for a zero 

reading of (P . - P ) , This means that the probe tip is aligned with the 
W 

flow direction in a horizontal plane. Then the five-way switching valves 
are set so that (P^ - Pg) is sensed by the pressure transducer. Finally, 
the reading of (P c - P w ) is similarly measured. 
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3 . 2 Calibra cions 

Two kinds of calibration are employed to reduce the raw data from 
the direct measurements. One is the calibration of the voltmeter, which 
determines a relationship between the voltmeter output and the magnitude 
of velocity. The other is the calibration of the five-hole pitot probe, 
which consists of two calibration characteristics: pitch angle 4 versus 

differential pressure ratio (P^ - Pg)/( P Q “ ? w ) » and velocity coefficient 

C - pV 2 /[2(P c - P w )] (3.1) 

versus pitch angle 5. 

Proper calibration is extremely critical to the accuracy of the ex- 
perimental work. The calibration of the five-hole pitot probe is based 
on the investigation, which concluded that the calibration coefficients 
are independent of Rep, which is based on probe tip diameter (5). The 
investigation determined that this conclusion is true for Rep > 1090, or 
a local velocity of 5.4 m/s. Low velocity might cause some calibration 
error. However, the investigation also shows that this calibration error 
affects the velocity measurements typically by less than 6 percent for 
Rep > 400, corresponding to a local velocity greater than 2.0 m/s. In 
this study, all calibrations are conducted at air velocity of 4.6 m/s. 

A schematic of the calibration facilities is given in Figure 8. The 
jet employed for calibrations is a free jet with a contoured nozzle which 
has a 3,5 cm diameter throat. Compressed air is used for the calibrating 
jet and its mass flow rate is controlled by a small pressure regular tor 
and a Fisher and Porter model 10A1735A rotameter. The five-hole pitot 
probe is rotated by the rotary table model BH-9 from Troyke Manufactur- 
ing Co., whose rotary vernier is readable to ± 0.5 minutes of arc. 


k „&jK3. . jb , 
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The calibration procedure of the voltmeter consists of reading the 

micro-manometer with a pitot static probe for a known air flow rate. 

3 3 

The volume flow rate is incremented by 0.5 ft /min [1 ft /min * 0.000472 
m'/s] over the range 0.5<Q<9.0. As shown in Figure 9, there exists a 
linear relationship between the velocity magnitude and the pressure read- 
ing. 

The first step to calibrate the five-hole pitot probe is to set the 
probe tip in the horizontal plane so that the yaw is aerodynamically 
nulled. This condition can be maintained by frequently checking the zero 
reading for (P tI - P„) during the whole calibration proceeding. The next 
step consists of reading the voltage output from the pressure transducer 
for differential pressures (P^ - Pg) and (P^ - P^) . These data are meas- 
ured at 5 degrees increments in 6 over the range -55 < 6 < 55. 

Figure 10 shows the calibration characteristics from which the pitch 
angle <5 and the velocity coefficient G are obtained, respectively. The 
calibration is conducted two times to be sure that the measurements are 
exact and that they remain the same. Even if the measurements are taken 
by different observers and at different times, the consistency is splen- 
did. Most of the discrepancies lie within 2-3 percent, and the maximum 
is within 10 percent. 

3.3 Data Reduction 

It is difficult to define in a word what turbulence is. According 
to Hinze f s (22) defninition, "Turbulent fluid motion is an irregular con- 
dition of flow in which the various quantities show a random variation 
with ti&Vr and space coordinates, so that statistically distinct values 
can be disconcerned." The vector velocity of turbulent flow consists of 
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a time-mean velocity and a fluctuating velocity component. As mentioned 

in the previous section, the differential pressure readings from the five- 

hole pitot probe are utilized directly to obtain the square of the vector 

velocity in the velocity coefficient. With the information from the five- 

hole pitot probe measurements alone, it is Impossible to express the 

“ 2 

square of the vector velocity V as the sum of the square of the magni- 

“2 “2 
tude of the time-mean velocity V and the time-mean fluctuation term V' , 

as given by 

V 2 » V 2 + V' 2 (3.2) 

2 

Obtaining the fluctuation term V' depends on other methods of meas- 

2 

urement. But the available data for the term V’ is not thought to be 
reliable. Furthermore, an established theory for the effect of turbul- 
ence in swirl flows on pressure probes is unavailable. So even though 
some discrepancy is to be expected, the reduced velocity is taken to be 
the time-mean velocity magnitude V which is written without the .overbar 
onward. In this study, all properties of the acting fluid are assumed 
to be those of ambient air. 

With the measurement data of the differential pressure ratio (P^ - 
P S )/(P C ~ P^j) > the corresponding pitch angle 6 is obtained ’with a cubic 
spline Interpolation technique from the calibration characteristics shown 
in Figure 10(a). The velocity coefficient C is determined from this 
value by using the corresponding calibration characteristic in Figure 
10(b). The magnitude of the velocity vector is calculated by 

v - [ f 0?c - V • C] 1/2 


(3.3) 
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The velocity components are easily calculated by elementary geometry 
using the measured yaw angle 3, the derived pitch angle <$, and the total 
velocity V. 

All data used in thi3 investigation were reduced with the FORTRAN 
data reduction computer program written by Rhode (5) . Full details of 
the data reduction procedure are given in Appendix C, including a program 
listing, sample input and output, and a user's guide. 



CHAPTER IV 


RESULTS 

Nonswirling and swirling nonreacting flows are investigated in an 
axisynunetric test section with expansion ratio D/d « 2, which may be 
equipped with a contraction nozzle of area ratio 2 and A. Velocity meas- 
urements are made with the five-hole pitot probe as described in Chapter 
lit. An analysis is made of the effects of various geometric parameters 
on the extent of the recirculation zones in the flowfield. These param- 
eters include side-wall expansion angle a * 90 and 45 degrees, swirl vane 
angle <j> = 0, 38, 45, 60, and 70 degrees, and contraction nozzle location 
L/D = 1 and 2 (if present) . The nozzle inlet velocities and Reynolds 
numbers employed in this study are high enough to ensure that the flow- 
fields are investigated under conditions independent of Reynolds number 
variation. All nozzle Inlet velocities and Reynolds numbers employed 
are listed in Table I. Flow characteristics are tabulated in terms of 
normalized u, v, and w velocity components, yaw angle 8 and pitch angle 
6 in Tables II through XII. Axial and swirl velocity profiles for all 
flows studied are shown in Figures 12 through 16 and 18 through 20. 

4.1 Effects of Swirl on Sudden Expansion 
Flows 

Swirling flows result from the application of a spiraling motion, 
with a swirl velocity component being Imparted to the flow via the use 
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of swirl vanes. A swirl number S characterizes the degree of swirl. The 
swirl vane angles employed in this study are 0 (swirler removed), 38, 45, 
60, and 70 degrees. These correspond to swirl numbers s. of 0, 0.52, 0.67, 
1.15, and 1.83, under the definition of swirl number introduced ia Section 
2.4. Figure 12 parts a through e show the axial and swirl velocity pro- 
files for $ » 0, 38, 45, 60, and 70 degrees, respectively, with side-wall 
expansion angle a * 90 degrees. 

The nonswirling flow investlaged is obtained with the swirler re- 
moved. Figure 12(a) shows an uniform axial velocity entering the test 
section. The corner recirculation zone extends to just beyond x/D = 2.0. 
The measurements for a corresponding flow taken with a stagnation tube 
and pitot tube by Chaturvedi (11). He found the reattachement point to 
be at x/D = 2.3 which is in good agreement with the present study. Moon 
and Rudinger (10) measured the reattahcment point with both theoretical 
and experimental methods in a similar circular test section with an expan- 
sion ratio D/d *» 1.43, which is different from the present study [with 
D/d = 2] . The result yielded a value of x/D = 1. 25 as a reattachment 
point, which corresponds to an attachment point approximately eight step- 
heights downstream. This is in good agreement with the present study. 

The velocity profiles for swirling flows shovm in Figure 12 parts 
b through e reveal that the flow entering through the swirl vanes is not 
uniform and has steep velocity gradients in the radial direction especial- 
ly at high swirl numbers. Furthermore, a considerable back flow around 
the hub is observed for = 70 degrees, as shown in Figure 12(e). For 
all values of swirl vane angle used in this study, the corner recircula- 
tion zone is not seen at x/D = 0.5 the closest axial location to the 
expansion block; instead, the maximum axial velocity is observed close 


22 


to the top wall at x/D « 0.5. The effects result from the strong centrif- 
ugal forces present in the incoming swirling flow. 

The central recirculation zone and precessing vortex core are now 
discussed. The precessing vortex core is defined as the region of high 
swirl, low axial velocity flow along the axis, which has a relatively 
constant small diameter. In flow visualization studied, it is seen to 
precess along the axis of the test section. The central recirculation 
zone is defined as the wide reverse flow region encountered near the in- 
let. Artistic impressions are given later in Figure 2 1 in which lines 
have been drawn connecting the radial positions of zero axial velocity. 

In case of the vortex core, that region is drawn along the zero axial 
velocity boundary in the downstream direction after the central recir- 
culation region* The size of the central recirculation zone increases 
with the increasing swirl vane angle until a certain value of swirl vane 
angle is reached (around 40 degrees) . Then its length begins to decrease 
under stronger swirl conditions. The results are tabulated in Table XIII. 
The core vortex was present at all values of swirl vane angle used in 
this investigation. In contrast to the central recirculation zone, the 
vortex core gets larger and larger continuously as the swirl vane angle 
increases. 

The swirl velocity peaks sharply around the edge of the expansion 
block before becoming more uniform farther downstream as shown in Figure 
12 parts b through e. A considerably nonuniform swirl velocity profile 
is observed at x/D - 0.5. Thereafter, relatively steady and uniform 
velocity profiles are seen, except for the region around the axis. The 
radial location where the maximum swirl velocity occurs goes up as the 
swirl vane angle increases. This trend is caused by the increase of 



23 


centrifugal effects. The swirl velocity along the axis is found to be 
zero as expected because of symmetry. 

4.2 Effects of Gradual Expansion on Flows 

Gradual expansion flows with a = 45 degrees were measured at two 
axial stations of x/D =0.5 and 1.0 for swirl vane angles of 0, 38, 45, 
60, and 70 degrees. Only the upstream flowfield needed to be thoroughly 
investigated in these cases, since inlet expansion effects affect this 
region the most, and their influence rapidly diminishes in the downstream 
direction (5, 13). Measurements were not taken at the inlet in this 
geometry because the presence of the expansion block interferes with 
probe positioning. 

The corresponding sequence of axial and swirl velocity ^profiles to 
a sudden expansion is given in Figure 13 parts a through e for a = 45 
degrees. Velocity profiles for a gradual expansion follow a similar 
trend to those for a sudden expansion. The major effect of a gradual 
inlet expansion is to encourage the air to flow along the side sloping 
wall, shorten the corner recirculation zone and accelerate axial velo- 
cities close to the top wall. Influence on the central recirculation 
zone for the swirl flow cases is minimal. 

4.3 Effects of Contraction Nozzle on Flows 

It is best to interpret the data obtained when the contraction block 
is inserted at different axial stations on the flowfields by comparing 
them to the data obtained without the contraction block. The measure- 
ments were taken at axial locations ranging from the inlet plane to the 
axial station just upstream of the station where the contraction block 
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was located. The effect of a contraction nozzle was investigated with 
the contraction block at L/D = 1 and 2 for a range of swirl strengths 
<f> = 0, 45, and 70 degrees with sudden expansion a = 90 degrees only. 
Figures 14, 15, and 16 present these velocity profiles with L/D = 1 and 
2 in parts a and b, respectively. 

Figure 14 shows that the contraction block generally affects the 
flowfield very little under nonswirling conditions. Furthermore, as the 
block is moved farther in the downstream direction, its effect on the 
flowfield decreases. The only effect that can be noticed from the velo- 
city profiles is decreased in length of the corner recirculation region 
due to the nozzle effect. 

Figure 15(a) shows the axial and swirl velocity profiles for the 
intermediate swirling flow of <j> = 45 degrees With a contraction block at 
L/D = 1. Significant positive axial velocities are observed along the 
center line. This situation is in striking contrast to that of the inter- 
mediate swirling flow without the contraction block, in which case the 
central recirculation zone spreads extensively along the centerline to 
x/D =1.5 with vortex core following it in the test section. Slightly 
larger swirl velocities are seen in Figure 15(a) with contraction block 
at L/D = 1 than in Figure 12(c) without contraction block. Velocity 
Profiels shown in Figure 15(b), for <j> *» 45 degrees with contraction block 
at L/D = 2, follow a trend similar to the corresponding flow without the 
contraction block at the first two axial stations of x/D = 0 and 0.5. 

But the axial velocities have quite different prof lies at two axial 
stations of x/D = 1.0 and 1.5 being affected by the proximity of the con- 
traction nozzle. At the axial station of x/D = 1.0, reverse flow is not 
observed at all, and considerable positive axial velocities are measured 
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near the axis. A rather uniform velocity profile is obtained at x/D * 
1.5. The effect of the contraction block on the swirl velocity profiles 
is negligible. The presence of a contraction block, in the intermediate 
swirling flowfield, results in the existence of positive axial velocities 
near the axis. 

For strong swirling flows of * 70 degrees, the axial and swirl ve- 
locity profiles are shown in Figure 16. A contraction block does not 
affect the flowfield much except at the axial station immediately up- 
stream of the block. 

In comparison with corresponding cases without downstream blockages, 
Figure 16(a) shows a wider central recirculation zone when L/D ■ 1 and 
Figure 16(b) shows a narrower precessing vortex core when L/D = 2. The 
swirl velocity profiles do not change significantly with or without 
presence of downstream blockages. 

In summary, a nozzle has most effect on the intermediate swirl case 
$ = 45 degrees. Its contraction effect in this case is strong enough to 
overwhelm the swirling recirculation region. 

However, a contraction has little effect on weakly swirling and 
strongly swirling flows, which are dominated by forward flow and centrif- 
ugal forces, respectively. 

Since the most drastic change in the velocity profile was observed 
at x/D =1.0 for intermediate swirling flow of <j> = 45 degrees with a con- 
traction block at L/D = 2, the measurements were repeated to ensure con- 
fidence in them. The new results, obtained at a slightly higher Reynolds 
number, confirm the previous results. Both sets of results are shown in 
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Figure 17, and seem to be very consistent. 

The effect of a stronger contraction nozzle was investigated for 
a range of swirl lengths * 0, 45 anJ 70 degrees with side-wall 
expansion angle a * 90 degrees. The contraction nozzle, of area ratio 
4 with 45 degree sloping upstream face, is located at 1/D * 1 and 2. 
Figures 18, 19, and 20 show these velocity profiles with L/D ■ 1 and 
2 in parts a and b, respectively. 

Figure 18 shows that the flowfield with a strong contraction nozzle 
changes very little as compared to the corresponding flowfield with a 
weak contraction nozzle shown in Figure 14. 

Figure 19 shows the axial and swirl velocity profiles for a swirl 
vane angle <p s 45 degrees. The presence of a strong contraction nozzle 
generates the high positive axial velocity near the axis. However, it 
deccelerates the axial velocity close to the top wall. The central re- 
circulation zone is smaller and located in an annular region, which 
excludes the axis. The swirl velocity profiles show narrower and 
stronger core swirl velocity gradients than previously. 

For swirl vane angle <J> = 70 degrees, the axial and swirl velocity 
profiles are given in Figure 20. It shows that the axial velocity near 
the axis is highly positive and the central recirculation region extends 
only to less than x/D = 1.0, much less than its corresponding case with 
the weak contraction block. At the axial station x/D = 1.0, forward 
flow occurs across the whole test section. This situation is in sharp 
contrast to that of a weak contraction nozzle, in which case the con- 
traction block affects the flowfiled very little under strong swirling 
conditions. Swirl velocity profiles given in Figure 20 show the 
narrower strong core swirl velocity gradient, as also observed in the 
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swirl vane angle (J> * 45 degree case. 

The above observances are summerized as follows. The effect of 
a weak contraction nozzle of area ratio 2 is confined to intermediate 
swirling flow case, while a strong contraction nozzle of area ratio 4 
affects both intermediate and strong swirling flow cases. Artistic 
impressions for a weak contraction nozzle and a strong contraction 
nozzle with emphasis on the recirculation zones are given in Figures 
20 and 23, for a range of swirl strengths with blockage located at 
L/D - 2. 


CHAPTER V 
CLOSURE 
5 . 1 Summary 

This study is five-hole pitot probe time-mean velocity measurements 
in confined swirling flows under low speed and nonreacting conditions. 
Three velocity components normalized with the inlet axial velocity are 
tabulated along with yaw and pitch angles. This provides an extensive 
data base for later verification of a theoretical simulation of the com- 
plex turbulent flowfield. 

The extent of the recirculation zone for sudden expansion flows with 
a * 90 degrees is characterized for swirl vane angles <j> = 0, 38, 45, 60, 
and 70 degrees. For a nonswirling flow, the corner recirculation zone 
extends to just beyond 2 with no central recirculation zone. The pres- 
ence of a swirler shortens the corner recirculation zone and generates a 
central recirculation zone followed by a precessing vortex core. The 
largest central recirculation zone was shown to be in the intermediate 
swirling flow case of $ = 38 degrees; the strongest precessing core vor- 
tex was found to be in the strongest swirling flow considered with <J> = 

70 degrees. 

Effects of a gradual inlet expansion of a » 45 degrees are investi- 
gated for the same swirl vane angles as employed for a sudden inlet ex- 
pansion. Velocity profiels for a gradual expansion follow a similar 
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trend to those for a sudden expansion. The only major effect of a grad- 
ual inlet expansion is to encourage the flow to remain close to the slop- 
ing side wall and shorten the extent of the corner recirculation zone in 
all cases investigated. 

Effects of the presence of a contraction nozzle at L/D *» 1 and 2 are 
Investigated for 4> * 0, 45, and 70 degrees. A contraction nozzle has 
little effect on weakly swirling and strongly swirling flows, which are 
dominated by forward flow and centrifugal forces, respectively. For the 
Intermediate swirl case of 4 1 “ 45 degrees, it encourages forward movement 
of otherwise slow-moving air and thereby shortens the central recircula- 
tion zone. A strong contraction nozzle of area ratio 4 has a more dramatic 
effect on the flowfields. 

5.2 Recommendations for Further Work 

The five-hole pitot probe technique is a useful cost-effective tool 
to investigate turbulent swirling recirculating confined flow. It en- 
ables time-mean axial, radial and swirl velocities to be deduced at any 
location in the flowfield. It is strongly recommended for use to obtain 
further flowfield measurements in the test facility at Oklahoma State 
University. It has, however, some inherent problems: 

1. Weak sensitivity to small velocities, which might permit large 
relative errors, 

2. Turbulence effects on the measurement of time-mean data cannot 
be easily accounted for. 

3. The measurement of normal and shear turbulent stresses is not 
possible. 

It is recommended that further work be done on the repeatability. 
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reliability, and accuracy of this technique. Results may be compared 
with hot-wire and laser doppler anemometers measurements in correspond- 
ing flow situations. 
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TABLE X 



NOZZLE 

INLET VELOCITIES AND 

REYNOLDS 

NUMBERS 



a = 

90° 

a - 

45° 

a = 90°, 

L/D = 1, 2 

<1* 

U in 

Re, 

a 

U in (rn/s) 

Re , 
a 

°in 

Re d 

0° 

15.7 

150,000 

15.5 

154,000 

22.3 

220,000 

38° 

10.5 

100,000 

10.6 

105,000 



45° 

12.6 

120,000 

14.9 

148,000 

13.5 

134,000 

60° 

8.84 

84,000 

9.58 

95,000 



70° 

5.57 

53,000 

6.25 

62,000 

6.8 

67,500 
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TABLE > 


CHARACTERISTICS OF FLOWF 
EFFECTS OF SWIRL VANE , 
EXPANSION ANG: 


Vane 

Angle 

(deg.) 

U. 

in 

(m/sec) 

Re d 

0 

15.7 

150,000 

38 

10.5 

100,000 

45 

12.6 

120,000 

60 

8.84 

84,000 

70 

5.57 

53,000 
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)S IN TERMS OF THE 
LE FOR SIDE-WALL 
a * 90° 


CTRZ 

(x/D) 

CRZ 

(x/D) 

W 

max 

(r/D) 


2.1 


1.6 


0.2 

1.5 


0.2 

1.3 


0.25 

1.2 


0.3 


l 

i. 

! 

I 

? 

I 
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TABLE XXV 

SAMPLE ALPHANUMERIC HEADINGS 


10000 

10040 

10050 

10060 

10070 

100S0 

10090 

10100 

10110 

10120 

10121 

10123 

10125 

10127 

10128 

10129 

10130 

10131 

10132 


COMPUTED MASS FLOW RATE 
COMPUTED MEAN AXIAL VELOCITY 
IJ VELOCITY 
V VELOCITY 
W VELOCITY 

TOTAL VELOCITY MAGNITUDE 
DIMENSIONLESS U VELOCITY 
DIMENSIONLESS V VELOCITY 
DIMENSIONLESS W VELOCITY 
DIMENSIONLESS TOTAL VELOCITY 

PROBE PITCH ANGLE < DEG ♦ ) 
PROBE YAW ANGLE (DEG.) 

P( NORTH) ~ p< SOUTH) (VOLTS) 

P( CENTER) - P ( WEST ) (VOLTS) 
MEAS. INLET MASS FLOW RATE 
MEAS. INLET AXIAL VELOCITY 
MEAS, INLET DYNAMIC PRESSURE 
MEAS. PITCH PRESSURE COEFF, 
AXIAL FLUX OF ANGULAR MOMENTUM 
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TABLE XV 

SAMPLE CALIBRATION DATA 


10133 

- 3*76 

10140 

- 2.81 

10150 

- 2*23 

10160 

- 1.72 

10170 

- 1.31 

10180 

- 1.00 

10190 

-.75 

10200 

-.59 

10210 

-.47 

10220 

-.35 

10230 

-.21 

10240 

-.063 

10250 

♦ 09 

10260 

.23 

10270 

.38 

10280 

.52 

10290 

.66 

10300 

.82 

10310 

1.00 

10320 

1 ♦ 23 

10330 

1 ♦ 49 

10340 

1.83 

10350 

2.24 

10360 

2.65 

10370 

3.4 


58 . 

2.42 

54 . 

1.86 

50 . 

1.59 

45 . 

1 . 35 

40 . 

1.19 

35 . 

1.11 

30 . 

1.02 

25 . 

♦ 95 

20 . 

.91 

15 ♦ 

.87 

10 . 

.85 

5 ♦ 

.88 

0 . 

.89 

•~5 ♦ 

.89 

—10 . 

.89 

- 15 . 

.89 

- 20 . 

.90 

- 25 . 

♦ 93 

” 30 . 

.96 

- 35 . 

1.04 

- 40 . 

1.10 

- 45 . 

1.24 

- 50 . 

1.40 

- 54 . 

1 . 59 

58 . 

2.04 
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TABLE XVI 

SAMPLE MEASUREMENT DATA 


10380 

0,0 

19 

,335 


10390 

0,0 

354,6 

-.007 

-.018 

10400 

, 3 

251,8 

♦ 003 

.021 

10410 

, 6 

258.0 

♦ 015 

.029 

10420 

,9 

262,6 

,016 

♦ 028 

10430 

1,2 

282,4 

.007 

.010 

10440 

1,5 

318.6 

-.004 

,205 

10450 

1,8 

324.6 

-.299 

2.450 

10460 

2,1 

329 . 0 

-.917 

4 . 308 

10470 

2.4 

326 . 6 

-1.010 

3.953 

10480 

2,7 

327.0 

-.865 

3 ♦ 255 

10490 

3,0 




10500 

3,3 




10510 

3,6 




10520 

3,9 




10530 

4,2 




10540 

4,5 




10550 

4,8 




10560 

5 , 1 




10570 

5,4 




10580 

5,88 

19 

,335 


‘ 10590 

0,0 

360 , 2 

,019 

,012 

10600 

.3 

300 . 0 

,036 

.007 

10610 

, 6 

274,4 

.060 

.051 

10620 

,9 

262 , 6 

♦ 081 

,130 

10630 

1,2 

256 , 0 

,094 

,213 

10640 

1.5 

o i=; o o 

,089 

.290 

10650 

1.8 

248,2 

. 088 

♦ 338 

10660 

2,1 

245,2 

.078 

.367 

10670 

2 , 4 

243.0 

♦ 067 

.352 

10680 

2 . 7 

241,6 

.059 

.324 

10690 

3.0 

243.0 

. 049 

,277 

10700 

3.3 

248.2 

,037 

,219 

10710 

3.6 

260 ♦ 0 

,035 

.138 

10720 

3.9 

282 ♦ 4 

, 038 

.093 

10730 

4,2 

308,6 

,046 

.156 

10740 

4,5 

325,4 

. 033 

♦ 338 

10750 

4,8 

332 . 0 

-,012 

.568 

10760 

5 ♦ 1 

334 . 2 

-.049 

,749 

10770 

5 .4 

331 . 0 

-.219 

. 727 

10780 

// 



END OF 

DATA 
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TABLE XVII 

SAMPLE AUXILIARY GEOMETRIC QUANTITY OUTPUT 


expansion angle*deg.i 


SWIRL VANE ANCLE* DEG. I 


9.009E+01 
4* 500E+0I 


inlet kadius*mi * 


combustor radius*mi = 


7.461E-02 
1 • 492E-01 


LAMINAR VI5C0SITY*KG/M/SECI 


1.600E-05 


1 • 16SE+00 


MEAS. INLET MASS FLO w RATE 


COMPUTED MASS FLOW RATE 


MEAS. INLET AXIAL VELOCITY 


DENSITY* XG/CU. Ml * 

““*•*”*-* — *-*-•» ***«•*••:?=•*—* 

1=1 2 

x = 0.0 0.14935 

2.76E-01 2.76E-01 

l = 1 2 

X = C.O 0.14935 

2.85E-01 2.91E-01 

*»«• 4- 4-*- *-*- *- *-# *•«.« 

1=1 2 

X = C.O 0.14935 

1. 35E+01 1.35E+01 

1=1 2 

X = 0.0 0.14935 

3 . 4 9E *-00 3. 57E4-00 

*-*-*-*. AXIAL FLUX OF ANGULAR MOMENTUM 

I = 1 2 

X = 0. 0 0.1 4935 

3.65E-C2 3.47E-02 


COMPUTED MEAN AXIAL VELOCITY 
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(b) . Simplified Test Section 

Figure 1. Typical Axisymmetric Gas Turbine Combustion Chamber 
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60 

30 


k Present data 
0 Reference data (5) 


-30 

-60 


- 1.5 - 1.0 - 0.5 0 0.5 1.0 

<p N - p s )/(p 0 - p w ) 

(a) Pitch Angle Calibration Characteristic 


k Present data 
0 Reference data (5) 


o6o° 


PITCH ANGLE 6 (DEG.) 

(b) Velocity Coefficient Calibration Characteristic 

Figure 10. Calibration Characteristics for Five-Hole 
pitot Probe 








(b) Swirl Velocity 


Figure 11. Flowfield Independence of Reynolds Number for 
Side-Wall Expansion. Angle a = 90° and 
Swirl Vane Angle <p = 45° at Axial Station 
x/D = 0.5. 






Figure 12. Velocity Profiles for Side-Wall Expansion Angle a = 90° without Contraction Block 






12 (Continued) 




Figure 12 (Continued) 



(a) u/u 
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Figure 13 (Continued) 








Figure 14 (Continued) 










Figure 16 (Continued) 
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. *, - ■»* 



o 


0.3 


U / Uo 


0.6 


(a) Axial Velocity 



Figure 17. Repeatability of Five-Hole Pitot Probe 

Measurement for Side-Wall Expansion Angle 
a = 90° and Swirl Vane Angle <j) = 45° with 
Contraction Block at L/D = 2,0 at x/D —1.0 




(a) u/u 
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Figure 18 (Continued) 
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x / D 


Figure 22. Artistic Impressions of Dividing Stream- 
lines with Contraction Dlock at L/D ; 
2.0 for Side-Wall Expansion Angle a ; 
90° and Swirl Vane Angle: (a) $ = 0° 

(b) $ = 45° (c) <j> = 70° . 
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0 1.0 2.0 3.0 


(b) 



(c) 


x/ D 

Figure 23. Artistic Impressions of Dividing Stream- 
lines with Contraction Block at L/D =2.0 
for Side-Wall Expansion Angle a = 90° and 
Swirl Vane Angle: (a) <j> = 0° (b) $ = 45° 

(c) cj) = 70°. 















APPENDIX C 


USER'S GUIDE TO FIVE-HOLE PITOT PROBE DATA 


REDUCTION COMPUTER PROGRAM 


A FORTRAN computer program to reduce five-hole pitot probe measure- 
ment data was written by Rhode (5) . The data reduction program consists 
of one main subprogram, two functions and five subroutines. A cubic 
spline interpolation technique is employed to interpolate calibration 
data. This user's guide emphasizes, the basic mathematical concepts 
together with analytical and numerical methods employed in this program. 
For user's convenience, sample input and output are given. The flow 
chart of this program is shown in Figure 24. 

1 . "MAIN Subprogram 

MAIN characterizes the particular conditions being investigated. 
Entire solution procedure is controlled by MAIN. So, it is the section 
of the program to which a user will devote most of his attention. It is 
divided into four chapters, each with a specific task, and a description 
of the individual chapters of MAIN now follows. 

CO Preliminaries 

Dimension and common are followed by input logical parameters, 
which activate, when specified as TRUE, certain special features of 
program. IWRITE writes solutions onto allocated disk storage and 
DIAGNS activates diagnostic write statements. Some geometric parameters 
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are specified as user input yalues. Input: read statement for alpha- 
nuraeric headings, is given. Initialization of all variables to be zero 
by calling subroutine INIT is' followed by two input read statements., 
which read calibration and raw measurement: data. 

Cl Data Reduction 

This is a main part to reduce raw measurement data. Data 
reduction consists of a cubic spline interpolation of calibration 
data and calculations of u, v and w velocity components. A pitch 
coefficient is calculated by (P^ t- Pg)/('P^ - P^) . With the calculated 
pitch coefficient, a pitch angle <$ is obtained from the corresponding 
calibration characteristic by calling function SPLINE. The correspond- 
ing velocity coefficient C for the resulting pitch angle 6 is also 
Obtained by calling function SPLINE once again. 

ID ID is a parameter to defihe caliration ranges. If pitch angle 
d or velocity coefficient € is out of the calibration range ,[-3.759 £ 

(P N - P s ) / (P c - P w ) £ 3,399’, |<S | £ 58], respectively, the correspond- 
ing value is, taken to be zero witli IDID = 0, 

With the above obtained information, three velocity components and 
magnitude of the total velocity vector are calculated from 


V « 1 1 (P c ~ P w ) * Cj 1/2 (A3.1) 

u - V cos. 6 cos 3 (A3, 2) 

v = V sin 6 (A3. 3) 

w = V cos 6 sin 3. (A3. 4) 
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C2 Auxiliary Calculations 


This part consist of calculations of geometric quantities and 
nondimensionizatlons of velocities. A mass flow rate and axial flux 
of angular momentum are calculated by a finite difference integration 
method as following; 


and 


mass flow rate = 2 u.Aa^ 

■i i i 


axial flux of angular momentum = 5 P u .w . r? A r 

i i i i x 

Mean axial velocity is calculated by 


(A3. 5) 


(A3. 6) 


mean axial velocity =' Eu^Aa^/A 

A nozzle inlet axial velocity is obtained from the conversion of 
the measured dynamic static pressure (inch H 2 0) at the nozzle throat 
into the corresponding velocity (m/ sec) . All velocities are nondimen- 
sionized by the nozzle inlet axial velocity. 

C3 Ouput 

All reduced data are printed out according to a standard format 
by calling subroutines WRITE and PRINT. 

2. INIT Subroutine 

In this subroutine, all variables are initialized to be zero 
even though some are set to obviously nonzero values. All variables 
in the subprogram MAIN communicate with those in the subroutine INIT 
through common statements. 
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3. SPLINE Function 


A cubic spline data fitting technique (25) is employed to inter- 
polate pitch angle <S and velocity coefficient C from the corresponding 
calibration characteristics. For a certain x value, which lies between 
the points (X^, Y^) and ^i+i> » a polynomial equation can be 

written by 

Y - a ± (X - X.) 3 + b^X - X i ) 2 + C t (X - X i ) + d ± (A3. 7) 

with application of the proper end point conditions, it reduces to 


h i-l S i~l + 2 (h i-l + h i' S i * h i S i+i 


b( 



Y - Y 
i i-1 


l i-l 


- } 


(A3. 8) 


h i = X i + 1 ' X i 


(A3. 9) 


During the foregoing procedures, the coefficients are also reduced 


as 


J i = (s i+i - 


b. = 


C. 

i 


S i /-2 

Y i+1 - Y i 
h i’ 


2h i S i + h i S i + l 


d. « Y. 
x x 


(A3. 10) 
(A3. 11) 

(A3. 12) 
(A3. 13) 


The equation (A3. 8) forms (n t- 2) equations for n unknowns. Two 
more equations needed to solve n Unknowns are obtained when the conditions 



>*» *v 
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pertaining to the end intervals are specified. Linear extrapolation 
techniquea is employed to obtain the end point conditions. The 
relations for end conditions of the whole curve are at the left end: 

h 2^ S 2 h l S 3 


h. 


(A3, 14) 


and the right end: 


S n = 


(h 0 + h n )S , - h _ S „ 
n~2 n-1 n-1 n-1 n-2 


(A3. 15) 


Finally, n equations Including two conditions are obtained for n 
unknowns. 

With the above mathematic concept in mind, one can follow the calc- 
ulation steps employed in this program. Two calibration and raw 
measurement data in MAIN communicate with those in SPLINE 
through parameters X, FX and XI. A Parameter A is an array for the 
coefficients of n simultaneous cubic polynomial equations. The (n-2) 
by (n-2) components of an array A are built by calling the subroutine 
ABUILD. The end point components A(1,J) and A(N,J) are built in 
the function SPLINE. The Gauss-Jordan elimination method is employed 
to solve (n-2) equations for (n-2) unknowns simultaneously. These 
solutions are translated back to the SPLINE by a parameter A. 

With the obtained solutions, the coefficient of a cubic poly- 
nomial equation are calculated by equations (A3. 10) through (A3, 13). 
Finally, an interpolation is achieved by substituting the known 
coefficient values into the equation (A3. 7). An interpolated value 
is translated back to the MAIN through a variable SPLINE. 


1 


4. H Function 


This function calculates the interval AX^ between the points 


v « i+ i- w- 


5 . ABUILD Subroutine 


This routine builds an (n-2) by (n-2) array for the coefficients 
of a cubi polynomial equation* This forms a tridiagonal matrix with 
a column vector of (n-*2) components. The notation employed here is as 
follows : 


a i,i-l " h l-l 


a i,i ■ ZO't-l + V 


a. 4 xi = b. 
x , 1+1 i 


,, Y i+i - Y i - Y i- Y t-i 1 
a i,n+l 6( h ± ' - 


6. GAUSS Subroutine 


This routine solves (n-2) equations simultaneously with a Gauss- 


Jordan elimination method. Solutions are stored in the a 


i,n+l’ 


7. PRINT and WRITE Subroutines 


The tasks of subroutines PRINT and WRITE are to pritn out the 
values of a 2-D array PHI* together with associated heading HEAD, the 


last parameter of the call list. The only difference between PRINT 
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and WRITE Is that WRITE is used to print out one line and PRINT is 
used to print out multiple lines at a time. The first two parameters 
of the call list, ISTART and JSTART, are usually 1 and 1, so that all 
internal and external values are printed to help diagnostics. 

8, Sample Input 

Alphanumeric headings, calibration data and raw measurement data 
are read by input read statements. Sample input data are given for the 
sample case of side-wall expansion angle a = 90 degrees, swirl vane 
angle <f> - 45 degrees and contraction nozzle located L/D = 1. 

(a) Alphanumeric headings: These alphanumeric headings are used 

to accompany the output descriptively. Nineteen lines of such headings 
are read in via the first read statement. The current alphanumeric 
headings are given in Table XIV. 

(b) Calibration Data: Calibrations are done for pitch angle and 

velocity coefficient characteristics. Calibration input data consist 

of pitch coefficient (P^ - P<J/(P C - P^) and velocity coefficient 
2 

pV / £ 2 (Pq - P^)] for 24 different pitch angles fi. Sample calibration 
input data are tabulated in Table XV. This table consists of pitch 
coefficient in the first column, corresponding pitch angle in the 
second and corresponding velocity coefficient in the third. 

(c) Measurement Data: Sample raw measurement data are given for 

two axial stations in Table XVI. Nineteen radial traverse measurement 
data are taken at each axial station. For each axial station, axial 
distance (inch), maximum number of radial positions and dynamic 
pressure (in E^O) are read as shown in lines 10380 and 10580 of Table 
XVI. Lines 10390 to 10570 are measurement data for the axial station 
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x/P * 0,0 and lines 10590 to 10770 are for the axial station x/D * 

0.5. The input measurement data consists of four column input data; 
the first is for radial measurement point (inch), the second for 
measured flow angle (360.0 - $ degree), the third for P N - Pg (volt) 
and the fourth for P G - (volt). 

9 , Sample Output 

Typical output consists of yaw angle, pitch angle and three 
velocity components. The corresponding output for the sample input 
data ;iust described is given in Table VI, as described earlier in 
Section 4,1 In this table, I and J denote the axial and radial station 
number, respectively. Also X and Y denote the normalized axial and 
radial station by the test section diameter, respectively. Yaw angle 
0 and pitch angle 6 are given in degrees. These velocity components 
are values normalized by the nozzle inlet axial velocity deduced 
from the pitot static pressure probe located upstream of the inlet 
to the test section, as described in Section 2.1, 

As auxiliary output, geometric specifications and quantities are 
printed out in front of the typical output. Geometric specifications 
include expansion angle (degree) . Swirl vane angle (degree) , inlet 
radius (m) and combustor radius (m) . On the other hand, flow 
quantities include mass flow rate (kg/s), nozzle inlet axial velocity 
(m/s), mean axial velocity (m/s) and axial flux of angular momentum 
(N.m). In each of these items, data relevant to each axial station are 
presented - these may be different with different axial stations and 
depend on the particular run conditions at the time of each tranverse. 
Output for the sample input data is given in Table XVII. 
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C SUBROUTINE MAIN 

C 

CA**********************************************T <" >■ '■’K****#*****## 

c 

c 

C A COMPUTER PROGRAM FOR BATA REDUCTION OF FIVE-HU..E PITOT 
C MEASUREMENTS IN TURBULENT* SWIRLING* RECIRCULATING* FLOW 

C IN COMBUSTOR GEOMETRIES 

C 

C VERSION OF JULY *1982 

C 

C H K YOON 

C MECHANICAL AND AEROSPACE ENGINEERING 

C OKLAHOMA STATE UNIVERSITY 

C STILLWATER* OK 74078 

C 
C 

CM************************************-*-**-*********************** 

CHAPTER 000000 PRELIMINARIES 00000000 
C 

DIMENSION HEDM < 9 > * HEDUMN ( 9 > *HEDNMS<9) *HEDCMW<9> * 

♦HEDU ( 9 > * HEDV ( 9 > * HEDW ( 9 > * HEDVT ( 9 > * HEDUST ( 9 ) * 

♦HEDVST ( 9 ) * HEDUST ( 9 ) * HEDVTS ( 9 > * HEDDEL < 9 > * HEDBET ( 9 ) * 

♦HEDMMF ( 9) *HEDMIV<9) *HEDMIPC9) *HEDMPP(9) *HEDAM(9) 

C 

COMMON 

♦/CALIB/CPITCH (24) *CDELTA(26) *CVELCF<26> 

t/MEASUR/RBETA (48*24)* RPNMPS (48*24)* RPCMPU (48*24) *NBATA ( 48 ) 


* * MAXJPT* RDNPRS ( 48 ) 

*/GE0M/X<48> *R<24> *XND(48> *RND(24> *BYPS(24> *DYNP(24> * 

* SNS(24) *NSTATN*XINCHS(4B> *RINCHS(24) 

*/CALC/VT0TAL(48*24> *U(48*24> *V(48*24) *W(48*24) * 

* VTSTAR(48*24> *USTAR(48*24) * VSTAR(48*24) * WSTAR(48*24) * 

* PICHCF (48*24)* VELCF (48*24)* DELTA (48*24)* BETA (48*24)* 
♦ANGMOM ( 48 ) * UMEAN ( 48 ) * MASS ( 48 ) * MASFLO ( 48 ) * UIN ( 48 ) 

REAL MASS* MASFLO 
LOGICAL IWRITE*DIAGNS 

C SET IURITE=.TRUE. FOR WRITING SOLN. ON DISK STORASE 

C SET DIAGNS=.TRIJE. TO ACTIVATE DIANOSTIC WRITE STATEMENT 

C — - NSTATN IS NO. OF AXIAL STATIONS FOR WHICH DATA IS 

C SUPPLIED 

C MAXJPT IS MAX. NO. OF RADIAL POSITIONS FOR ALL AXIAL 

C STATIONS 

C CPITCH IS CALIBRATION PITCH COEFF. 

C CDELTA IS CALIBRATION PITCH ANGLE(DEG. ) 

C CVELCF IS CALIBRATION VELOCITY COEFF. 

C — NCAL IS NO. OF CALIBRATION POINTS 

C XINCHS IS AXIAL PQSITION( IN. ) OF EACH PARTICULAR RADIAL 

C TRAVERSE 

C — — NDATA IS NO. OF RADIAL LOCATIONS FOR WHICH DATA IS 

C SUPPLIED FOR EACH PARTICULAR RADIAL TRAVERSE 

C RDNPRS IS MEASURED DYNAMIC PRESSURE ENTERING THE 

C SUIRLER 

C RBETA IS MEASURED FLOW ANGLE (DEG.) WHERE YAW 

C ANGLE BETA « 360.0 - RBETA 

C RDNPRS IS MEASURED VOLTS FOR PNORTH - PSOUTH 

C DIFF. PRESSURE 

C — RPCMPW IS MEASURED VOLTS FOR PCENTER - PWEST 

C DIFF. PRESSURE 

C PICHCF IS REDUCED VALUE FOR PITCH COEFF. 

c DELTA IS REDUCED VALUE FOR PITCH ANGLE (DEG. ) 

C VELCF IS REDUCED VALUE FOR VELOCITY COEFF. 

C BETA IS REDUCED VALUE FOR PROBE YAW ANGLE (DEG* ) 

c VTOTAL IS TOTAL VECTOR VELOCITY MAGNITUDE(M/SEC) 

c U IS AXIAL VELOCITY <M/SEC) 

C — V IS RADIAL VELOCITY (M/SEC) 

C W IS SWIRL VELOCITY (M/SEC) 

C — — VTSTAR IS DIMENSIONLESS TOTAL VELOCITY 

C — ' USTAR IS DIMENSIONLESS AXIAL VELOCITY 

C VSTAR IS DIMENSIONLESS RADIAL VELOCITY 

C WSTAR IS DIMENSIONLESS SWIRL VELOCITY 
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C ALU PRIMARY USER INPUTS ARE LOCATED HERE 

I WRITE’. TRUE. 

DIAONS-, FALSE. 

ALPHA-90 . 

PH1=4S, 

VISCOS=l . 8E-5 
NSTATN=2 
MAX JPT=19 
PATM=75.1 
TATM=20. 

IT-- 4 8 
JT=24 

RLARGE=T\ 75/ <2. 0i«39. 37) 

RSMALL=RLAR'GE/2 . 0 

READ ( 5 f 205 ) HEDM f HEDUMN t HEDU . HEDV f HEDW f 


* HEDVT , HEDUST f HEDVST f HEDUST f HEDVTS , HEDDEL t HEDBET r 
I H EON MS . HEDCMU . HEDMMF , HEDMI V f HEDMIP f HEDMPP f HEDAM 
205 FORMAT (9A4) 

C INITIALIZE VARIABLES TO ZERO 

CALL INIT 

C READ FIVE-HOLE PITOT CALIBRATION DATA 


NCAL=25 

DO 10 I=1fNCAL 

READ ( 5 f 210 ) CPITCH < I ) f CDELTA ( I > f CVELCF < I ) 

10 CONTINUE 
210 FORMAT (3F10 .5 ) 

IF (DIAGNS) WRITE < 6 f 400 ) (CPITCH ( I) f 1=1 >25) 

IF (DIAGNS) WRITE ( 6 f 400 > (CDELTA( I ) fI=1 f25) 

IFCDIAGNS) WRITE(6f400> ( CVELCFC I ) f 1=1 f 25) 

400 F0RMAT(///f1Xf13(F8.4f1X> f//f5Xf 12(F8.4) ) 

C READ RAW MEASURED DATA TO BE REDUCED 

DO 30 I=1fNSTATN 

READ! 5 f 230 ) XINCHS ( I > f NDATA ( I ) f RDNPRS ( I ) 

JPTS=NDATA( I ) 

DO 20 J=1fJPTS 

READ <5 f 220) RINCHS ( J ) f RBETA ( I f J > f RPNMPS < I v J ) f RPCMPW ( I .■ J ) 
20 CONTINUE 
30 CONTINUE 

C CONVERT X AND R FROM INCHES TO METERS 

DO 35 I=1fNSTATN 
X< I )=XINCHS( I)*0.0254 
JPTS=NBATA< I > 

DO 32 J=1fJPTS 
R < J) -RINCHS ( J ) )X0 ,0254 
32 CONTINUE 
35 CONTINUE 
220 FORMAT (4F10.5) 

230 FORMAT ( IF 10 .5f1I, IOfIFIO.5) 

IF ( DIAGNS ) WRITE < <f> f 470 ) ( NDATA ( I ) f 1=1 f NSTATN ) n 
IF (DIAGNS) WRITE ( 6 f 450 ) ( X ( I ) f 1 = 1 f NSTATN ) 

IF (DIAGNS) WRITE (6f 500) < R ( J ) f J=1 f JPTS ) 

DO 37 1=1 f NSTATN 

IF (DIAGNS) WRITE ( 6 f 500 ) < RBETA ( I f J ) f J=1 f JPTS ) 

IF ( DIAGNS ) WRITE ( 6 f 500 ) (RPNMPS ( I f J) f J=1 f JPTS ) 

IF (DIAGNS) WRITE ( 6 f 500 ) ( RPCMPW ( I f J) f J=1 f JPTS ) 

37 CONTINUE 

450 FORMAT (///f40Xf1(F8.4f1X) ) 

500 F0RMAT(///f20Xf10(F8.4) ) 

C 

CHAPTER 11111 DATA REDUCTION 111111 
C 

470 FORMAT ( /// f 40X f 1 ( 18 f IX) ) 


C CALC PICHCF AND INTERPOLATE FOR DELTA FROM 

C PITOT CALIBRATION CURVE 


IDID-0 

DO 50 1=1 f NSTATN 
JPTS=NDAT A ( I ) 

DO 40 J=-1fJFTS 

IF ( ( RPCMPW ( I f J ) .EO. 0.0) .AND. ( RPNMPS ( I f J )-. EO . 0. 
♦ GO TO 38 
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PICHCF < t » J > =RPNMP8 (I f J > / ( RPCMPW (I f J > + 1 . E-A > 

IF( (PICHCF < f f J) *GT .3*399) .OR. (PICHCFt I , J) .LT.-3.7S9) > 

* GO TO AU 

IFdDID .EO. 0) OKLTAd f J)=SPLINE(CPITCHf 
I CDEL r A fNCALf PICHCF < I f J > ) 

IFdDID .GT. 0) DELTACIf J)=SP(CPITCH f CDELTAr 

* NCALfPICHCF(IfJ) > 

IDID-1 

GO TO 40 
38 CONTINUE 

I.iF.LTA(If J)"0,0 
WRITE(6 f850) IfJ 

850 F0RMAT(20Xf 'PICHCF IS OUT OF RANGE OF CALIBRATION AT 1= 

#' fI3f' ANI) J»' r 13) 

40 CONTINUE 
SO CONTINUE 

C INTERPOLATE FOR VELCF FROM PITOT CALIBRATION DATA 

IDID-O 

DO 80 I=1fNSTATN 
JPTS-NDATAd ) 

DO >0 J=-1f JPTS 

IF ( ( RPCMPW d f J) .EG. 0.0) .AND. <RPNMPS(IfJ> .EO. 0.0>) 

I GO TO A5 

IF(<ABS(OELTA<IrJ>) > .GT. 58.0) GO TO <45 
IFdDID .EG. 0) UELCFdf J)=SPLINE(CDELTA f 
iCVELCF f NCAL f DELTA (I f J > ) 

IF < IDID .GT. 0) VELCFdf J)==SP(CDELTAfCVELCFf 
*NCALf DELTA (If J) ) 

IDID-1 
GO TO 70 
65 CONTINUE 

VELCF(IfJ)=0.0 
WRITE(6f8?0> IfJ 

890 FORMAT (20X» 'DELTA IS OUT OF RANGE OF CALIBRATION DATA 
♦AT 1=5' f 13 r ' AND J='fI3) 

70 CONTINUE 
80 CONTINUE 

DO 85 I=1fNSTATN 

IF CDIAGNS) WRITE(6f 500) ( PICHCFd f J ) f J=1 r JPTS) 

IF(DIAGNS) WRITE<4f 500> (DELTAdrJ) fJ=1fJPTS> 

IF(DIAGNS) WRITE (6 f 500) (VELCFd r J) f J»1 t JPTS) 

85 CONTINUE 

CALC MAGNITUDE OF TOTAL MEAN VELOCITY VECTOR AND 

Uf y, & U COMPONENTS 

RHO=PATM* d.36 . 0/0 . 1 02 ) / ( 287 . 0* < TATM+273 . 0 ) ) 

PI=3 . 1-5159 
DO 100 1=1 f NSTAT N 
JPTS” NDATAd) 
no 90 J -If JPTS 
BET A( I f J) -360,-RBCT A < I f J) 

IF((RPCMPW(IfJ).E 0,0.0) .AND. (RPNMPSd f J) .EG. 0.0) ) 

I BETA(IfJ)=0.0 

VT 01 ALdFJ)=SaRT(AEiS(2.0/RH0*VELCFdF J)*RPCMPW(If J) * 133.9) > 
UdFJ)=VTOTALdFJ) * COS ( DELTA ( I f J)*PI/180.0) * 

* COS < BETA ( I f J)*PI/180 .0 ) 

Vd f J)=V'liiTAL (IfJ) * SIN ( DELTA (I r J)*PI/180. 0 ) 
U(IfJ)=OTOTAL(IfJ> * COS(DELTA(If J)*PI/180.0) * 

t SIN(BETA(IfJ>*PI/100.0) 

90 CONTINUE 
100 CONTINUE 

IF (DIAGNS) WRITE(6»500> (VTOTALd r J) fJ=1fJPTS) 

IF ( DIAGNS ) WRITE < 6 f 500 > CU (I f J) f vJ=l f JPTS > 

IF (DIAGNS) WRITE ( 6 f 500 > ( V ( I> J) fJ* IfJPTS) 

IF (DIAGNS) WRITE ( 6 r 500 > ( W(I f J) fJ -If JPTS) 

CHAPTER 222222 AUXILIARY CALCULATIONS 22222 
C 

DO 130 I=1fNSTATN 

C CALC GEOMETRIC QUANTITIES 

JPTS=NDATA( I ) 
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JPTSM1=JPTS-1 

DYPS(1)=0.0 

DYNP < JPT3 ) =2 . 0* < RL ARGE-R ( JPTS ) ) 

DO 110 J*1fJPTSM 1 
DYNP<J>=R<J+i>-R< J> 

DYPS<J+1)«DYNP<J> 

110 CONTINUE 

DO 115 J=1fJPTS 

SNS < J > =0 , 5* < DYNP ( J > +D YPS ( J ) ) 

115 CONTINUE 

IF (DIAGNS) WRITE(6f500) <DYNP( J) f J=1f JPTS) 

IF ( DIAGNS) WRITE(6fS 00) <8NS( J) r J=>1, JPTS) 

FL0U=0,0 
UMOMmO.O 
DO 120 J=1f JPTS 

WMQM=WMOM+RHO*R ( J > **2*SNS < J > *U < I f J ) *W < 1 f J ) 

FLQW=FL0W+2 . 5»<PI*R ( J ) *SNS ( J > *IJ (I f J ) 

IF (DIAGNS) WRITE(6f?00) Jf FLOW f WHOM fRHO 
900 FORMAT <///f 1'5f4<F12.4f2X) ) 

120 CONTINUE 

ANGMOM<I)=WMOM 

MASSCI)=RHO*FLOW 

UMEAN < I > =MASS < I ) / ( RHO*PI *RLARGE**2 > 

130 CONTINUE 

IF(DIAGNS) WRITE(6f 450) (UM C AN(I) f 1=1 fNSTATN) 

IF (DIAGNS) WRITE<6f450) CMASS(I) f 1=1 fNSTATN ) 

C NONOIMENSIONALIZE VELOCITIES 

DO 150 1=1 fNSTATN 

XND < I ) =X < I > / < 2 . 0*RLARGE ) 

JPTS=NnATA(I> 

UIN (I ) = ( SORT < 2 . 0/RH0*RDNPRS (I>*249.08) )#(6» 312/5 . 94 > ** 
MASFLO ( I > =2 . 0*PI*RH0#UIN ( I > *RSMALL**2/2 . 0 
DO 140 J=1fJPTS 

VTSTAR<IfJ)=VTOTAL<IfJ)/UIN<I) 

USTAR< I f J)=U( I f J)/UIN ( I ) 

vs tap < i f j > =v < i f j ) /uin a ) 

WSTAR<IfJ)=WCIfJ)/UIN<I) 

140 CONTINUE 
150 CONTINUE 

DO 160 J=1fMAXJPT 
RND< J)=R< J)/C2.0*RLARGE> 

160 CONTINUE 
C 

CHAPTER 33333 OUTPUT 333333 3 

C 

l'FT.NOT. IWRITE) GO TO 165 
WRITE(ll) X 
WRITE (11) R 
WRITE <11 > U 
WRITE(ll) V 
WRITEdl) W 
165 CONTINUE 

WRITE<6f311 ) 

WRITE(6f325) ALPHA 
WRITE ( 6 f 330 ) PHI 
WRITE<6r335> RSMALL 
WRITF.<6 f340) RLARGE 
WRITE(6f355) VISCOS 
WRITE(6f360> RHO 

c 

CALL WRITE < 1 f 1 f NSTATN • 1 f IT t JT f X f R f MASFLO f KEDMMF > 

CALL WRITE ( 1 f 1 f NSTAI N f 1 f IT r JT fXf R f MASS f HEDM ) 

CALL WRITEUfIfNSTATNfIfITf JTfXfRfUINfHEDMIV) 

CALL WRITE<1f1fNSTATNf1fITfJTfXfRfUMEANfHEDUMN> 

CALL WRITEC 1 f 1 fNSTATN fIfITfJTfXfRf ANGMOMf HEDAM) 

CALL PRINT < 1 f 1 f NSTATN f MAX JPT f IT f JTfX f R f U f HEDU > 

CALL PRINT ( 1 f 1 fNSTATN f MAX JPT f ITf JT fXfRf VfHEDV) 

CALL PRINT <1 , 1 , NSTATN f MAX JPT f IT f JT f X f R f W f HEDW ) 

CALL PRI NT ( 1 f 1 f NSTATN f MAX JPT f IT f JT F X f R r DELTA f HEDDEL ) 
CALL PRINT < 1 ? 1 f NSTATN f MAX JPT f IT f JT f X f Rf DETA f HEDPET ) 
CALL PRINT (1 f1 fNSTATN f MAX JPT f ITf JT fXfRf VTOTALf HLDVT > 
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CALL PRINT < 1 f 1 f NSTATN . MAX JPT f IT f JT f XND r RND f USTAR f HEBUST > 
CALL PR INT ( 1 f 1 f NSTATN f MAX JPT r I T f JT f XND f RNC! f VSTAR f HEBVST ) 
CALL PR INT < 1 f 1 r NSTATN f MAX JPT t IT f JT r XND f RND t WSTAR f HEDWST ) 
CALL PRINT< 1 f 1 f NSTATN » MAX JPT f ITf JT fXND fRND f VTSTAR fHEDVTS ) 
CALL PR INT < 1 f 1 f NSTATN f MAX JPT f IT f JT f XINCHS f RINCHS f RPNMPS 

* rHEDNMS) 

CALL PR INT ( 1 f 1 , NSTATN f MAX JPT f IT f JTf XINCHS f RINCHS f RPCMPU 

* fHEDCMU) 

CALL WRITE <1 f 1 f NSTATN f 1 f IT f JT f XINCHS fRINCHS f RDNPRS f HEDMIP ) 
CALL PRINT ( 1 f 1 f NSTATN f MAX JPT * IT f JTf XINCHS fRINCHS fPICHCF 

* fHEBMPP) 

STOP 

c FORMAT STATEMENTS 

311 FORMAT < 1H1 fT37f 'AXISYMMETRICf ISOTHERMALf GT COMBUSTOR 

* FLOWFIELB MEASUREMENTS' f//fT53f 'USING A FIVE-HOLE 

* PITOT PROBE') 

325 FORMAT < //// fT40 r 'EXPANSION ANGLE < DEG. ) =' fT77f1PE13.3> 

330 FORMAT <//fT40f 'SWIRL VANE ANGLE (BEG.) =' fT77f 1PE13.3) 

335 FORMAT <//fT40f 'INLET RADIUS(M) =' fT77f1PE13.3> 

340 FORMAT (//*T40r 'COMBUSTOR RADIUS(M) = ' fT77f 1PE13.3) 

355 FORMAT <//fT40f 'LAMINAR VISCOSITY<KG/M/SEC > =' fT77f 1PE13.3) 
3A0 FORMAT (//fT40f 'DENSITY (KG/CU. M) =' fT77f 1PE13.3f////> 

END 

C 

SUBROUTINE INIT 

C****************)t*********************************************** 

C 

COMMON 

IVMEASUR/RBETA < 48 f 24 ) » RPNMPS <48 f 24 ) f RPCMPW ( 48 f 24 ) r NDATA < 48) 

* f MAX JPT f RDNPRS (48) 

*/GE0M/X<48> fR<24 ) » XND ( 48 ) fRND ( 24) f DYPS< 24 ) r DYNP(24 ) t 

* SNS<24) f NSTATN fXINCHS< 48) f RINCHS < 24) 
l/CALC/VTOTAL<40r24)fU<48r24) fV<48f24) fW(48f24) f 

I VTSTAR < 48 f 24) t USTAR <48,24 ) f VSTAR < 48 f 24) » WSTAR <48.24 ) t 

t PICHCF< 4 Bf24 ) f VELCF < 48f24 > » DELTA < 48 f 24 ) f BETA < 48 f 24 ) f 

* ANGM0M<48> fUMEAN<48> fMASS<48) fMASFL0<48) fUIN<48) 

REAL MASSrMASFLO 


DO 20 I=«l, NSTATN 
MASFLO< I )=0.0 
MASS ( I )=0 . 0 
ANGMOM( I )~0 . 0 
UMEAN < I ) “0 .0 
UIN<I>»0.0 
DO 10 J«=1 , MAXJPT 
VTOTAL< I f J)=0.0 
U<If J)~ 0.0 
V<If J)=0.0 
W(If J)=0.0 
VTSTARCIf J)=0.0 
USTAR < I » J)=0.0 
VSTAR (If J>=0.0 
WSTAR< I f J) =0 . 0 
RBETA < I f J)=0 . 0 
BETA < I f J)=0.0 
RPNMPS (I fJ)=0.0 
RPCMPWdr J)*0.0 
PICHCF < I f J)=0 .0 
VELCF < I f J)=0.0 
DELTA < I f J) =0.0 
.10 CONTINUE 
20 CONTINUE 
RETURN 
END 


FUNCTION SPLINECXf FXf Nf XI) 

C*****************m*****#**********************»*********»***** 

CUBIC SPLINE CURVE FITTING IN 2 DIMENSIONAL DATA PLANE 
INPUT VALUES S 

Xr FX DATA ARRAYS t ONE DIMENSIONAL r X IN INCREASING ORDER 
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C N NUMBFR OF DATA POINTS IN X? MAX 2& 

C XI POINT OP INTEREST ? WHERE F(Xi> IS TO BL FOUND 

C 

C RETURN VALUE i 

C SPLINE OR SP » F CXI ) 

C THIS ROUTINE ACTIVATES ROUTINE ABUILD? Hr AND GAUSS. 

C FOR INTERPOLATION OF A LARGE NUMBER OF DATA POINTS? FUNCTION 

C SPLINE MAY BE CALLED ONLY ONCE ? AND SUBSEQUENT CALLS MAY US 

C ENTRY POINT SP. 

DIMENSION X < 1 > r FX(l)? A(2A?27) 

C CONSTRUCT SPLINE MATRIX 

NF'1=N+1 



DO 10 1=1 r 

N 


DO 10 J~1f 

NF‘l 

10 

A< It J)=0« 



NM1-N-1 



DO 20 1*2 ? 

NH1 

20 

CALL ABUILDCX? FX? A? 


A<1f1)«H<Xf 

2) 


A < 1 f2)~-H<X 

f1)-H<Xf2) 


A ( 1 f 3 ) =H ( X f 

1) 


NM2-N-2 



A ( N ? NM2 ) =H C X r NM i ) 

A < N ? NM1 ) ~-H ( X ? NM2 > -H ( X ? NM1 ) 

A(N?N)-H(X?NM2) 

C FIND SECOND DERIVATIVES 

CALL GAUSSCA? N? NP1 ) 

ENTRY SPCX? FX? N? XI) 

C FIND FCX1) 

DO 40 1=1? NM1 

IP1=IF1 

ircxi . EO . '.(D) GO TO SO 

IFCXl .LT, XCI) .AND. XI , GT. XCIPl)) GO TO 41 
IF <Xl .GT. XCI) .AND. XI .LT. XCIPl) ) GO TO 41 

40 CONTINUE 

1FCXJ .EO. XCN) ) GO TO AO 
WHITt.CA? 42) XI 

42 FORMAT (' Xl=', G14.7? ' OUT OF INTERPOLATION RANGE? 

I RETURNED VALUE=0 ' ) 

SP~0. 

SPLINE- 0. 

STOP 

41 CONTINUE 
IP1-I+1 
HI”H(X? I > 

HX“X1-X C I > 

AI = (AC IP1 ?NP1 )-AC I?NP1) )/ < A . #HI) 

HI a A C I ? NP1 ) /2 . 

O'l CFX(IPl>-FXCI>>/HI-C2.*HI*ACI?NPl)+HI*A(IPi?NPl>>/4. 

DI = FXCI) 

SPLINE=AI*HX**3+BI#HX**2+CI*HX+DI 

SP--SPLINE 

RETURN 

C 

50 CONTINUE 

SPLINE=FXCI> 

SP=SPLINE 

RETURN 

C 

AO CONTINUE 

SPLINE=FX(N> 

SP=SPLINE 

RETURN 

END 

C 

FUNCTION HCX? I ) 

c***M*mmmm*******m*w***********x%***************m******** 

C CALCULATE DELTA X WHICH IS USUALLY CALLED AS H. 
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DIMENSION X<1> 

IP1=I+1 

H=X<IP1>-X<I> 

RETURN 

END 

C 

SUBROUTINE ABUXLDCXf Ff Af Nf I) 

C ******!m******************************************************* 
C CONSTRUCT SPLINE MATRIX FOR FINDING 2ND DERIVATIVES. 

C**H0tt*5K#*^10tt^***)*t**********)tt3*C*)*C*)»C)tt****)»t)*t)tt)tt*»)tt*)*t**J*(******3»C**J»tm)*C*)»t<t 

DIMENSION X<1>f FC 1 > f A(26f27> 

IH1=I-1 

IP1=I+1 

NP1=N+1 

HI=H<XfI> 

HIM1«H<XfIM1> 

A < I f IM 1 )=HIMl 
A<IfI>=2.*<HIM1+HI> 

A<IfIP1)=HI , , 

A(IfNP 1>=< <F<IPi)-F<I>)/HI - <F(I>-F<IM1))/HIM1 )* 6 • 

RETURN 

END 

C 

SUBROUTINE GAUSS<Af Nf NP1> 

P***^*****^****^***********#******************************* 
C GAUSS-JORDAN ELIMINATION 

i. ******im*********¥*##**»^)M*****#*#****##**#***********##**#*** 

DIMENSION A<2Af27> 

NM1 = N-l 

DO 3 L=1f NM1 

LP1=L+1 

DO 3 I=LP1fN 

CONST=AUfU/A(LfL> 

DO 3 J-Lf NPi 

3 A(IfJ)=A(If J)-CONST*A<Lf J> 

DO 6 1=1 f NM1 

IP1=I+1 

DO 6 L=IP1fN 

CONST=A<IfL)/A<LfL> 

DO 6 J=I f NPI 

& A<I.J>»A(Ir J)-CONST#A<L fJ> 

DO 10 1=1 f N 

A(IfNP1)®A(IfNP1>/A<IfI> 

10 A<IfI)»1. 

RETURN 

END 

C 

SUBROUTINE PRINT< ISTARTf JSTARTfNIfNJfITf JTfXf YfPHI fHEAD) 
CA#*****#***#**************#*******#*************)m************* 
c 

DIMENSION PHKITfJT)fX(IT) fY< JT> rHEAD<9 > fST0RE<48> 

ISKIP=1 

JSKIP=1 

URITE<Af110)HEAD 
ISTA= I START-12 

100 CONTINUE 
ISTA=ISTA+12 
IEND=ISTA+11 

IF(NI . LT • I END) IEND=NI 

WRITE(6f 111 > ( I f I=ISTAf IENDf I SKIP) 

URITE<6f114) CX( I > f I=ISTAf IENDf ISKIP) 

WRITE<6f112) 

DO 101 JJ=JSTARTfNJfJSKIP 

j=jstart+nj-jj 

DO 120 I=ISTAfIEND 
A=PHI<If J) 

IF(ABS(A) .LT.l.E-20) A=0.0 
120 STORE! I )=A 

101 WRITECAf 113) Jf Y< J) f <STORE(I)fI=ISTAf IENDf ISKIP) 

IF< IEND . LT . NI >G0 TO 100 

RETURN 
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110 F0RMAT<1H0f17<2H*-)f7Xf9A4f7Xf17(2H-*>> 

111 FORMATUHOf 13H I = fI2f11I9> 

112 FORMAT (BHO J Y> 

1U FORMAT < 13 fOPF8 .5f IX f 1P12E9 .2) 

114 FORMAT <11H X = fF8.5f 11F9.5) 

END 

SUBROUTINE URITE< I START f JSTART fNIfNJfITfJT fXf YfPHI fHEAD) 
fl****#*##*##**####*********#*****##******##**##********#******* 

DIMENSION PHI < IT) f X < I T ) f Y ( JT > f HEAD < 9 > fSTORE<48) 

ISKIP 1 *! 

JSKIP»1 

WRITE < 6 t 110) HEAD 
I8TA'=ISTART-12 

100 CONTINUE 
ISTA=ISTA+12 
IEND-ISTA+11 

IF « NI . LT . 1 END ) IEND=NI 
WRITE<6f111> (IfIsISTA^IENDfISKIP) 

WRITE< 6f 1 14 ) < X< I > f I=ISTAf IENDf I SKIP) 

DO 101 JJ= JSTART fNJf JSKIP 
,1-JSTART+NJ-JJ 
HO 120 IsISTA fIEND 
A- PHI < I > 

ir<Al'S(A).LT.l.E-20) A=0 . 0 
120 STORE* I) -A 

101 WRITE ( 6? 1 13 ) (STORE (I) fI e ISTAfIENDfISKIP) 

IF < IENDfLT . NI >80 '10 100 

RETURN 

110 FORMAT < 1H0 f 17 (2H*-) f 7Xf9A4 f7Xf 17 <2H-*> ) 

111 FORMAT* lH0fl3H I » fI2f11I9> 

113 FORMAT < 12Xf 1P12E9* 2) 

114 FORMAT <11H X = fF8.5f 11F9.5) 

C 

END 
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SUMMARY 

The present study clarifies swirler performance characteristics as an 
aid to understanding and computer modeling of gas turbine combustor flowfields. 
The swirler under study is annular with 10 adjustable vanes and a hub to 
swirler diameter ratio of 0.25. Measurements of time-mean axial, radial, 
and tangential velocities at the swirler exit plane for a range of blade 
swirl angles <p from 0 to 70 degrees show clearly the effects of centrifugal 
forces, recirculation zones, andbladestall on the combustor inlet velocity 
profiles. 

INTRODUCTION 

The problem of optimizing gas turbine combustion chamber design is 

complex, due to the many conflicting design requirements. In an effort to 

reduce development time and cost, computer simulation techniques are being 

used to study the effect of geometry changes on the flowfield inside combustor 
1 2 

designs, ’ Accuracy of predictions from a computer model is strongly dependent 
on the inlet boundary conditions used, which are primarily determined 
by the swirler and its performance at different vane angle settings. In 
previous studies,^ the velocity boundary conditions at the inlet to the model 
combustor were approximated by idealized flat profiles for axial and swirl 


velocity, with radial velocity assumed to be zero. However, measurements taken 
close to swirl er outlets show that the profiles produced are quite nonuniform 
as found elsewhere 3 and at Oklahoma State University. 

In the present study, the objective was to make velocity measurements as 
close as possible to the swirl er exit, so as to define better the performance 
characteristics of the swirler. A range of swirl-blade angles $ from 0 (no 
swirler present) to 70 degrees is being considered, with 0, 45 and 70 degrees 
being hiqhlighted. Previous velocity measurements at the inlet to the test 
section on the same test facility^ were obtained with expansion blocks attached 
to the exit of the swirl pack, thus being 3.2 cm downstream of the swirler exit 
plane. In the present study, these blocks were removed so as to allow 
measurements directly at the swirler exit plane and give a better indication 
of the actual performance of the swirler. 

EQUIPMENT AND PROCEDURE 

The swirler is annular with a hub one-fourth the inside diameter of the 
swirler housing, with hub and housing diameters of 3.75 and 15.0 pm, 
respectively. The 10 vanes (or blades) are attached to shafts which pass through 
the housing wall and allow individual adjustment of each blade's angle. The 
blades are wedge-shaped to give a constant pitch-to-chord ratio of approximately 

5 

one, which should give good turning efficiency. 

The time-mean velocity components are measured with a five-hole pitot 
probe, which allows determination of the magnitude and direction of the mean 

velocity vector simultaneously. The probe is mounted in a traversing 
mechanism which allows it to be translated vertically and rotated about the 
yaw axis. Tubing from the probe's five pressure taps is routed through 
selector valves to a differential pressure transducer, and the resulting 
pressure difference values are read directly from an integrating digital 
voltmeter. The pressure data are reduced by a computer program to yield 


nondimensionalized u, v, and w velocity components, which are then plotted 
in the form of profiles. Details appear in Ref. 4. 

RESULTS 

Velocity profiles for a number of flowfields have been compiled for both 
radial and aximuthal traverses, and some sample results are presented below: 

1 , Radial Profiles 

Figure 1 shows radial profiles of axial and swirl velocity running 
from the swirl er centerline to twice the swirl er exit radius 
[r/D =0.5 where D is the test section diameter used in associated 
studies] for three values of blade angle <p : zero (no swirler), 

45, and 70 degrees. The nonswirling case has a nearly-flat axial 
velocity profile, as expected. All velocities shown are normalized 
with respect to the mean exit axial velocity, measured independently. 
However, for the 45 degree case, both profiles vary from zero at the 
hub to a maximum at or near the rim of the swirler in an almost linear 
fashion. This is due to centrifugal forces on the rotating air, 
resulting in much of the flow leaving near the outer boundary. 

Note the low velocities for r/D < 0.1 approximately; they are 
probably an effect of the central recirculation zone extending upstream 
enough to obstruct the inlet flow. 

The 70 degree profile shows strong reverse flow in the range in the 
central region showing that the central recirculation zone now extends 
upstream past the exit plane. The maximum velocity at the outer edge 
is approximately three times the inlet velocity and the velocity 
gradients across the profile are quite large. The positive axial 
velocity adjacent to the hub is probably the result of the slight 
clearance between the blades and the hub, allowing a thin layer of 
nonswirling air to pass directly through the swirler. 


2. Azimuthal Profiles 


Figure 2 shows the azimuthal or 0-variation of axial, radial 
and swirl velocities for the case of <J> = 45 degrees. Similar 
results occur at other swirl angles. The profiles span an angle of 
48 degrees, more than the 36 degrees between successive blades. 

The measurements were taken at a constant radial distance of 
r/D = 0.20. 

The profiles all show significant variation with azimuthal position, 
which can be attributed to several causes. Among them are blade 
stall from using flat blades at a high angle of attack, not enough 
blade overlap (pitch/chord ratio too high), and wakes from blunt 
trailing edges. The radial component is seen to be of the same order 
of magnitude as the axial and swirl components, and also shows the 
greatest variation. 

CONCLUSION S 

Assumptions of flat axial and swirl profiles with radial velocity 
equal to zero are progressively less realistic as the swirler blade 
angle <j> increases. At moderate swirl = 45 degrees, linearly increasing 
profiles of u and w with radius are appropriate with v nonzero. At 
strong swirl <f> = 70 degrees, even more spiked profiles are appropriate 
with most of the flow leaving the swirler near its outer edge. Nonaxi symmetry 
prevails in all swirl cases investigated. 
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(a). <p = 0° (No Swirler) 
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(b) <f> = 45° 
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(c) 4> = 70° 


Fig. 1 Radial Profiles of Axial and Swirl Velocity at Swirler Exit Plane 
for Swirl Cases: (a) 4> = 0°, (b) <J> = 4B°* and (c) 4> = 70°. 
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(b) Radial Velocity 
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(c) Swirl Velocity 


Fig. 2 Azimuthal Profiles of (a) Axial Velocity, (b) Radial 

Velocity, and (c) Swirl Velocity at Swirl er Exit Plane 
for Swirl Case d> = 45°. 
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SUMMARY 

The present contribution provides numerical predictions of turbulent 
swirling recirculating confined flows using various inlet velocity starting 
conditions for the case of swirl vane angle equal to 45 degrees. The validity 
of flowfield predictions resulting from the choice of inlet profiles is 
assessed by comparing the predicted velocity profiles with corresponding 
experimental velocity profiles. Results demonstrate that realistic 
predictions are forthcoming only from the inclusion of realistic axial, radial 
and swirl velocity profiles as inlet conditions, and that considerable errors 
occur if unrealistic idealized inlet conditions are used. Predictions are 
then exhibited for a range of swirl strengths [<j> = 0, 38, 60 and 70 degrees] 
using measured inlet axial, radial and swirl velocity profiles in each case. 
The ensuing flowfields are characterized via velocity profiles and streamline 
patterns, and illustrate the large-scale effects of inlet swirl on flowfields. 
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1, INTRODUCTION 

The governing differential equations of swirling recirculating flows are 
elliptic, and solutions depend strongly on the boundary conditions applied 
around the flow domain. It is important to define adequately the boundary 
conditions, especially the inlet velocity profiles. Predictions given in Ref. 
1 have been obtained by using flat inlet axial and swirl velocity profiles in 
which the swirl velocity magnitude is deduced from the swirl vane angle <j>. In 
Refs. 2 and 3 the initial values of the components of mean axial and swirl 
velocity were taken from experimental data, while the radial component of mean 
velocity was assumed negligible. On the other hand, Ref. 4 assumes the swirl 
velocity profile to be that of solid body rotation. 

The present contribution provides numerical predictions of the velocity 
field for a swirling confined jet flow through a vane swirler, expanding into 
a larger chamber with diameter ratio D/d = 2. A swirl vane angle of 45 
degrees is selected to evaluate, via comparison with time-mean flow 
measurements, the predictions using various types of inlet axial, radial and 
swirl velocity profiles. This illsutrates their effect on the flowfield, and 
dramatizes the need to use realistic inlet boundary conditions in predictive 
studies and turbulence modeling improvements. Also, predictions of the' 
velocity field for a complete range of swirl vane angles (0, 38, 60 and 70 
degrees) using the best measured inlet velocity profiles are exhibited. They 
illustrate the large-scale effects of inlet swirl on flowfields. 
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2. TYPES OF INLET BOUNDARY CONDITIONS CONSIDERED 

In all cases considered, the kinetic energy of turbulence k and its 
dissipation rate are specified as in general accepted ways. 6 The total 
time-mean velocity magnitude has been used in specifying k. Four possible 
specifications of the inlet velocities are considered: 

1. Flat inlet axial and swirl velocities with radial velocity zero are 
assumed. That is, both u and w are constant valued; 
u Q = constant 
w Q - u 0 tan <{> 

where is the best estimate of the constant flow turn angle, 
slightly less than the swirl vane angle. [When blades = 45 deg. ^ 
has been found experimentally that <J>fi ow = 42 deg.]. 

2. As Case 1., except that the inlet swirl velocity profile is assumed 
to be that of solid body rotation: 
u Q = constant 

w o = w m,o * r /R 

where w mj0 is the maximum orifice value of w which occurs at the 
outer edge r * R of the inlet. The value of w m)0 is so chosen as to 
maintain the swirl number 5 the above as Case 1 for -simi larity 
purposes as is now described. Defining 

6 = w m,o 7u o 

it can be shown 6 that S and G are related via 
S = (G/2 ) / [1 - (G/2) 2 ]. 

For a given flat profile flow angle <P = <!>fi ow , it can also be 
shown 6 ’ 7 that approximately: 


S = 2/3 tan <j>. 


Thus the equivalent swirl number jet with different inlet profile 
types can be deduced. 

Measured inlet axial and swirl velocities are used® with radial 
velocity assumed to be zero. 

Measured inlet axial, radial and swirl velocity values are used, 
taking data from recent five-hole pitot-probe data in close vicinity 
of the swirler exit.® The measured data are linearly interpolated to 
fit into the computer code as shown in Table 1. 
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3, RESULTS AND DISCUSSION 

Predicted results are obtained via the STARPIC computer program, ^ 
developed at Oklahoma State University as an advanced version of the Imperial 
College, London, TEACH computer code. 1 -® The present version includes the 
calculation of swirl velocity, and incorporates its effect on wall 
functions. Currently the standard k-e turbulence model is employed, with 
minor changes as in Ref. 1. Sudden and gradual expansions can be accommodated 
easily. It has already been used' 1 - for tentative predictions using the 
profiles of Case 1 described in Section 2. 

Predictions are first discussed which deal with the flow through inlet 
swirl guide vanes of angle 4> = 45 degrees, which impart a flow angle of 
approximately $ = 42 degrees to the flow, although this varies a little with 
radial location. Various types of inlet profile assumptions are considered 
and their similarities and differences are noted. Predictions are then 
presented for flowfields with four other swirl strengths using measured inlet 
flow velocity profiles for axial, radial and swirl velocity in each case. To 
aid convergence, the solution for each value of <P [0, 38, 45, 60 and 70 deg.] 
is used for the initial starting values for the next higher value of (j>. [This 
includes the predictions for the specific case of <j> = 45 deg. discussed at 
length in Section 3.1]. The convergence criteria in each case are all 
normalized residual source sums to be less than 0.004 or a maximum of 200 
iterations, with results in Table 2 and Figs. 1 through 19 being so generated. 
3.1 Effect of Inlet Velocity Profiles 

Figure 1 shows the predicted velocity profiles at various downstream 
axial stations using flat inlet profiles for axial and swirl velocities with 
inlet radial velocity equal to zero. This is Case 1 of Section 2. For 
purpose of comparison with the experimental data, the averaged flow angle 


value is used instead of <J> = 45 degrees. The cited value in Ref. 8 is 42 
degrees. A central recirculation zone is seen to be more than one chamber 
diameter in length, and approximately one half of the chamber diameter in 
diameter, encouraged by centrifugal force. No corner recirculation zone is 
observed beyond x/D = 0.5. 

Figure 2 displays the predicted velocity profiles using flat inlet axial 
velocity profile with the inlet radial velocity is equal to zero. The inlet 
swirl velocity profile is assumed to be that of solid body rotation for flow 
angle of 42 deg. This is Case 2 of Section 2. A flow recirculation is 
observed at the corner of the combustor through x/D = 0.5 with a very small 
recirculation region at the center. 

Flowfield prediction using measured inlet axial and swirl velocity 
profiles 7 with the assumption that inlet radial velocity is equal to zero are 
displayed in Fig. 3. Two recirculation zones are predicted, one on the corner 
and the other being a larger zone located on the axis. The central zone 
extends downstream up to approximately 1.5 combustor diameters, with maximum 
core diameter of 0.4 times the chamber diameter. Notice that this is somewhat 
longer and narrower, with higher reverse velocities, than found in the 
previous case. The corner zone extends to less than one diameter in length. 

Predictions made by using the most realistic measured inlet axial, radial 
and swirl velocity profiles taken from Ref, 8 are shown in Fig. 4. The corner 
recirculation zone terminates before x/D = 0.5 and a wider central 
recirculation zone than previously calculated is observed. The flowfield 
structure is illustrated via streamline plots, which are computer calculated 
and drawn, as shown in Fig. 5. A three-dimensional representation of the 
predicted axial velocity is given in Fig. 6, so as to aid in visualizing flow 


zones . 
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Figure 7 displays the measured axial and swirl velocity profiles using a 
five-hole pitot probe for the corresponding inlet flow swirl case.^ A larger 
central recirculation zone which extends downstream more than two chamber 
diameters is observed. The flow fills the larger chamber rapidly because' of 
centrifugal effects, and no corner recirculation zone is seen in the velocity 
profiles shown, although flow visualization studies reveal a corner zone of 
length approximately 0.4 chamber diameters. From comparison of the predicted - 
velocity profiles with the measured velocity profiles, it can be seen that 
flowfield predictions following the most realistic inlet velocity 
specifications give best agreement with the data, see Figs. 4 and 7. 

3.2 Effect of Swirl Vane Angle on Flowfield Patterns 

The velocity field predictions for swirl vane angles of 0, 38, 60 and 70 
degrees using the obtained more realistic inlet velocity profile are displayed 
in Figs. 8, 9, 10, and 11. The streamline patterns are computed and plotted 
as shown in Figs. 12, 13, 14, and 15 and the flowfield construction is shown 
in Figs. 16, 17, 18, and 19. .The following results can be drawn from these 
figures: 

1. The higher the swirl strength, the higher the centrifugal force, hence 
the larger the central recirculation zone core diameter; 

2. The stronger the swirl strength, the longer the central zone extends 
downstream; 

3. The corner recirculation zone size reduces sharply as swirl strength 
increases 

Recirculation zone lengths are- dramatically affected as the degree of 
swirl increases. Discrepancies seen on comparing predictions with 
corresponding flowfield data obtained with a five-hole pitot probe^ indicate 
typically longer corner and center zones predicted than measured. This may be 

6"3 
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because of poor probe sensitivity iri turbulent low velocity regions, and/cr 
poor turbulence model performance in these regions. Only further detailed 
hot-wire and/or laser doppter anemometer measurements and turbulence model 
development will resolve the inconsistencies. Despite this, the results are 
in general agreement with the discussion of Ref. 7. 

3.3 Convergence Rates 

Table 2(a) displays the residual source sums (RSS) after 200 iterations 
between successive values of tj>. It can be seen that as the swirl vane angle 
increases, the RSS increases. This is expected because of the complexity of 
the flow at the higher swirl strengths. In Part (b) of the table, the effect 
of inlet velocity profile assumption on the degree of convergence for the <|> = 
45 deg. case is illustrated via the RSS values. It can be seen that as the 
degree of nonun if ormity in inlet profiles increases, the degree of convergence 
reduces, with higher RSS values after a given number of iterations. This is 
generally to be expected. 

Table 3 displays the progress toward convergence with more stringent 
final convergence criteria [alt normalized residual source sums to be less 
than 0-004 or a maximum of 500 iterations]. Part (a) shows normalized 
residual source sums for u, v, w and mass; Part (b) shows nondimensionalized 
velocities u/u 0 , v/u Q and w/u 0 at a monitored field location (1,0) = (8,8). 
Both parts of the table contain the sequence of five degrees of swirl. It is 
found that the convergence rate for <Ji s 0 deg. case is good enough at 200 
iterations, The degree of convergence in all the cases can be clearly seen. 
However, Figs. 20 and 21 give u and w profile plots for the cases <J> = 45 and 
70 degrees, respectively, after complete convergence, requiring 315 and 500 
iterations, respectively. Comparison of these figures with Figs. 4 and 11, 
which were obtained after 200 iterations between successive swirl cases, 
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reveals that general field results were good enough after 200 iterations, 

I 

although the fine detail and turbulence data may clearly require quite 
stringent convergence criteria. 

Figures 22 and 23 display the prediction of the flowfield streamlines 
pattern for the cases $ = 45 and 70 deg., respectively, obtained after 
complete convergence, requiring 315 and 500 iterations, respectively. Again, 
they are comparable with Figs. 5 and 15 which were obtained with a lower 
number of iterations between successive swirl cases investigated. 

3 . 4 Effect of Inlet Turbulence Energy 

Table 4 displays the predicted axial and swirl velocities and turbulence 
parameters for swirl vane angle <l> = 45 deg. at the axial location x/D = 0.5. 
Only the weak degree of convergence was used for those results requiring 200 
iterations in each case. Three types of specification of inlet kinetic energy 
of turbulence k were considered: 

1. Inlet k is assumed to be a constant obtained from the average exit 

axial velocity: 

k = 0.03 a y ,o 2 

2. Inlet k is allowed to vary with radius and is a function of local 
inlet time-mean axial velocity at each radial location 

k - 0,03 u 0 2 

3. Inlet k is allowed to vary with radius and is a function of local 
inlet toal time-mean velocity at each radial location 

k = 0.03 (u 0 2 + v 0 2 + w 0 2 ) 

The factor 0.03 in these equations is the recommended constant given inlet 
relative turbulence intensity of 0.245. Specimen predictions using these 
three inlet k assumptions with realistic inlet velocity values are given in 
Table 4, Parts a, b, and c, respectively. They illustrate that there is 



little change in velocity predictions with changes in the inlet k 
specification. 
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4, CONCLUSIONS 

A numerical prediction using the standard two-equation k- e turbulence 
model and different inlet flow boundary condition assumption has been applied 
to a confined swirling flow. Inlet flow boundary conditions have been 
demonstrated to be extremely important in simulating a flowfield via numerical 
calculations. Predictions with either flat inlet profiles, solid body 
rotation or zero radial velocity are inappropriate. Predictions are given for 
a full range of swirl strengths using measured inlet axial, radial and swirl 
velocity profiles in each case, and predicted velocity profiles and streamline 
plots illustrate the large-scale effects of inlet swirl on flowfields. 
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Fig. 16 Axial velocity representation for nonswirling flow. 

Fig. 17 Axial velocity representation for swirl vane angle <f> = 38°. 

Fig. 13 Axial velocity representation for swirl vane angle <j> = 60°. 

Fig. 19 Axial velocity representation for swirl vane angle <p = 70°. 

Fig. 20 Predicted velocity profile" for swirl vane angle <J> = 45° using 

measured inlet u, v and w profiles and high degree of convergence. 

Fig. 21 Predicted velocity profiles for swirl vane angle <f> = 70° using 

measured inlet u, v and w profiles and high degree of convergence. 

Fig. 22 Predicted streamlines for swirl vane angle 4> = 45° using measured 
inlet u, v and w profiles and high degree of convergence. 

Fig. 23 Predicted streamlines for swirl vane angle <j> - 70° using measured 
inlet u, v and w profiles, and high degree of convergence. 
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fig. 5 Predicted streamlines for swirl vane angle <j> = 45° using 
measured inlet, u, v and w profiles. 
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Fig. 6 Axial velocity representation for swirl vane angle $ = 45° 







Predicted velocity profiles for swirl vane angle $ « 38° using 
measured inlet u, v and w profiles. 










AXIAL POSITION X/D 


Fig. 13 Predicted streamlines for swirl vane angle <)> = 38° using 
measured inlet u, v, and w profiles. 

















Fig. 17 Axial velocity representation for swirl vane angle $ = 38° 





Fig. 13 Axial velocity representation for swirl vane angle 


ORIGINAL Bit B 

OF POOR QUALITY 


€ 

10 

» 

0 



Fig. 19 Axial velocity representation for swirl vane angle 4> = 70 
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Fig. 22 Predicted streamlines for swirl vane angle <{> = 45° using measured 
inlet u, v and*w profiles and high degree of convergence. 
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Fig. 23 Predicted streamlines for swirl vane angle <J> = 70° using measured 
inlet u, v and w profiles, and high degree of convergence. 
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Table 1. Fitting Experimental Data into the Computer 
Code for Swirl Vane Angle <t> - 45 deg. 


Wall Expansion Angle ot = 90 deg. 

Average Inlet Velocity u Q = 13.0 m/s 

(a) Axial Velocity 

Experimental Data: J r/D u/u Q 


11 

0.25532 

0.000 

10 

0.22979 

1.610 

9 

0.20425 

1.570 

8 

0.17872 

1.450 

7 

0.15319 

1 .030 

6 

0.12766 

0.618 

5 

0.10213 

0.327 

4 

0.07660 

0.332 

3 

0.05106 

0.000 

2 

0.02553 

0.000 


J 

r/D 

u/u 0 

u (m/s) 

8 

.2425 

0.8085 

10.51 

7 

.215 

1.5866 

20.63 

6 

.175 

1 .3890 

18.05 

5 

.13 

0.6560 

8.52 

4 

.085 

0.3304 

4.29 

3 

.045 

0.000 

0.00 

2 

.0125 

0.000 

0.00 


Interpolated Data: 


Table 1 (continued) 

(b) Radial Velocity 
Experimental Data: 


Interpolated Data: 


J 

R/D 

v/u 0 

10 

0.22979 

0.5450 

9 

0.20425 

0.4890 

8 

0.17872 

0.5050 

7 

0.15319 

0.5090 

6 

0.12766 

0.5690 

5 

0.10213 

0.3590 

4 

0,07660 

0.3420 

3 

0,05106 

0 .0000 

2 

0.02553 

0.0000 


J 

R/D 

v/u 0 

v(m/s) 

8 

0.22875 

0.5427 

7.05 

7 

0.19500 

0.4948 

6.43 • 

6 

0.15250 

0.5106 

6.64 

5 

0.10750 

0.4032 

5.24 

4 

0.06500 

0.1867 

2.43 

3 

0 .02875 . 

0.0000 

0.00 

2 

0.00000 

0.0000 

0.00 


Table 1 (continued) 
(c) Swirl Velocity 


Experimental Data: J r/D 

11 0.25532 

10 0.22979 

9 0.20425 

8 0.17872 

7 0.15319 

6 0.12766 

5 0.10213 

4 0.07660 

3 0.05106 

2 0.02553 


Interpolated Data: J r/D 

8 0.2425 

7 0.2150 

6 0.1750. 

5 0.1300 

4 0.0850 

3 0.0450 

2 0.0125 


. t 


w/u 0 

0.000 
1.410 
1.450 
1.320 
0.862 
0.580 
0.508 
0.186 
0.000 
0.000 

E. 

i 

| 

jj 

w/u 0 w (m/s) j 

} 

■■ - - ■ ■ |i 

0.7080 9.20 ! 

1.4334 18.63 

1.2533 16.29 

0.6058 7.88 

0.2919 3.80 

0.0000 0.00 

0.0000 0.00 





Table 2. Absolute Residual Source Sums After 200 Iterations 


(a) For different swirl vane angles, with measured u, v and w 


Swirl Vane Angl e <f> 

u 

V 

w 

Mass 

0 

0.00560 

0 .00096 

0.00500 

0.00009 

38 

0 .05060 

0.02340 

0.13190 

0 .04280 

45 

0 .02040 

0 .02200 

0.02580 

0.00720 

60 

0.03607 

0.02960 

0.04286 

0.01410 

70 

0.05322 

0 .05032 

0.09408 

0.04622 

(b) For different inlet profile assumptions 

for <f> - 45 degrees 

Inlet Type 

u 

V 

w 

Mass 

Flat u and w with v zero 

0 .00634 

0.0048 

0.0039 

0.0012 

Flat u with w solid body 
rotation and v zero 

0.02540 

0.0270 

0.0313 

0.02700 

Measured u and w with 
v zero 

0.03509 

0.0393 

0.0426 

0.01055 

Measured u, v and w 

0 .02040 

0 .0220 

0.0258 

0.00720 



w 


Mass 


0.00096 


0.005 


0.0001 


0.0234 

0.0029 

0.0042 

0.0009 


0.1319 

0.0218 

0.0053 

0.0039 


0.0428 

0.0073 

0.0027 

0 . 002 ' 


0.0190 

0.0040 

0.0037 


0.0183 

0.0048 

0.0039 


0.0058 

0.0012 

0.0005 


0.0296 

0.0131 

0.0065 

0.0037 


0.0482 

0.0249 

0.0126 

0.0065 


0.0159 

0.0085 

0.0041 

0.00156 


0.0352 

0.0177 

0.0079 

0.0044 


0.0694 

0.0317 

0.0151 

0.0071 


0.0230 

0.0177 

0.0079 

0.0042 






































Table 3. (Continued) 

(b) Nondimensional Velocities at Field Location (I* J) s (8, 8) 


0 

<f> 

Iteration 

Number 

U/U 0 

< 

< 

O 

W/U 0 

0 

200 

0.705 

- - - : 

0.0126 

0.0 


200 

0.680 

0.035 

mm 

38 

300 

0.645 

0.027 

Bn 

400 

0.663 

0.030 



430 

0.665 

0.031 

in 


200 

0.199 

0.056 

mm 

45 

300 

0.170 

0.058 



315 

0.169 ■ 

0.058 

mam 


200 

0.022 

0.074 

— 

60 

300 

0.044 

0.071 

mmm 

400 

0.055 

0.070 ! 

BUM 


500 

0.060 

0.067 

0.646 


200 

0.185 

1 

0.063 

1.064 

70 

300 

0.143 

0.066; 

1.024 

400 

0.116 

0.068 

1.010 


500 

0.105 

0.069 1 

1.00 




Table 4. Effect of Inlet Turbulence Parameters on Flowfield Predictions 


45° and at x 


(a) k = 0.03 u g (constant) 


.365 
.529 
0.003 


.49 

9050.0 

26.5 

.606 

10939.0 

32.4 

.635 

21025.0 

76.2 


(b) k = 0.03 u Q 2 (varies with radius) 


0.4 

0.377 

0.487 

7960.0 

22.8 

0.3 

0.535 

0.606 

8832.0 

25.0 

0.2 

0.026 

0.652 

19462.0 

69.8 


(c) k = 0.03 (u 2 + v 2 + w 2 ) (varies with radius) 


0.377 

0.542 

0.026 


7460.0 

8090.0 

19537.0 
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Abstract 


Numerical finite difference predictions are 
made of recirculation zones behind bluff-body flame 
stabilizers, showing quantitatively the effects of 
forebody geometry, blockage ratio, lateral position 
of the blockage and inlet swirl on the central re- 
circulation zone. A simple transient Navier-Stokes 
solution algorithm and laminar flow simulation are 
used with * free slip* and *no slip 1 wall boundary 
conditions, thus illustrating how a basic approach 
may be used to solve a sophisticated fluid dynamic 
problem. 

Introduction 

The Phenomenon 

Flames in combustion chambers are generally 
stabilized by means of recirculation zones. 1 10 
The recirculation zone acts as a continuous heat 
source by transporting hot combustion products up- 
stream towards the burner and thereby supplying the 
energy required for ignition* In addition to the 
stability of the flame, the size and shape of the 
recirculation zone affects the volumetric heat re- 
lease rate and the combustion efficiency. 

Previous experimental work shows that the re- 
circulation zone can be set up in the wake of a 
bluff-body such as a stabilizer disk or cone. 
Alternatively it can be generated by the action of 
swirl and then, at sufficiently high degrees of 
swirl, the presence of a central bluff-body is not 
required. In a divergent nozzle, where the flow 
remains attached to the wall, the increased radial 
component of the flow results in an increase in 
the size of the recirculation eddy. The combined 
effect of blockage, swirl and divergent nozzle 
considerably enhances the process. This paper is 
concerned with nonreacting flows; a schematic of 
the flow system appears in Fig. 1. 

The Problem 

To aid the design of combustors with complex 
geometries, one wishes to characterize this recir- 
culation zone theoretically, both qualitatively 
and quantitatively. The properties of this zone 
(length, breadth, reverse mass flow rate, and maxi- 
mum reverse flow velocity) provide a key to the 
residence time distribution and time-temperature 
history of gases in the system, r\\d bear a strong 
relation to efficiency of, rad ollution cmergir. 
froffly t h e syst em . 
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ing, Member AIAA. 

**Associate Professor, Dept, of Mechanical 
Engineering 

***Professor , Dept, of Mechanical Engineering 
****Associate Professor, School of Mechanical and 
Aerospace Engineering, Associate Fellow AIAA. 

Copyright (e) American Institute of Aeronautics and 
Astronautics. Inc., I9S2. All rights reserved. 



Fig. 1 Schematic of the axisymmetric problem 
with disk blockage 

Mathematically the recirculating flow requires 
the solution of the full Navier-Stokes equations 
and a transient time-marching technique is used 
here. The present investigation omits combustion 
and treats the turbulence of the flow field rather 
crudely: the flow is supposed to behave as a lami- 

nar flow with a constant artificially high value 
ascirbed to the viscosity. Even so, it is interest- 
ing to see what level of accuracy results from this 
quite reasonable assumption. 

The Present Contribution 

Consideration is given here to a primitive, 
pressure- velocity variable, finite difference code 
which has been developed to predict 2-D axisymmetric 
swirling transient flows. The technique used here 
incorporates the following: 

(i) A finite difference procedure is used in 
which the dependent variables are the 
velocity components and pressure, formu- 
lated in axisymmetric cylindrical polar 
coordinates. 

(ii) The pressure is deduced from the continu- 
ity equation and the latest velocity 
field, using a guess-and-correct proce- 
dure for the latter. 

(iii) The procedure incorporates displaced 
grids for the two velocity components, 
which are placed between the nodes Where 
pressure is stored. 

The code is designed for persons with little 
or no experience in numerical fluid dynamics with 
the purpose of demonstrating that many useful and 
difficult problems can be attacked without resorting 
to large, complicated computer programs and useful 
qualitative computer experimentation may be per- 
formed. As specific complexities can be added 
easily, the program provides a basis for developing 




many new numerical capabilities* The present work 
has grown out of previous ideas about solving 2-D 
flows using transient Los Alamos 11 and steady 
state Imperial College 14 15 methods. 

Solution Technique 
The General Idea 


A laminar flow simulation with 'free slip' 
or 'no slip' wall boundary conditions is used 
in simplified code, based directly on the SOLA 
algorithm 13 , suitably extended to include swirl 
effects. It provides an alternative to more 
sophisticated codes that are_also applicable to 
the problems of this type. 14 *” 15 It is assumed 
that the flow is 2-D isothermal, incompressible 
and axis ymme trie, and chat material properties 
are constant. Under these assunptions , the 
governing equations can be expressed as follows: 
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(uvr) + (vw) + - V(V*w-^j) 


where 


v j * - + |4 + - I 4 

y 3x 3y* x 3x 


Here u, v and w are velocity components in the 
x(E r), y(S symmetry axis) and 0 directions, respec- 
tively, in accordance with the basic algorithm. 13 

A general grid system covering the domain of 
interest appears in Fig. 2, illustrating how bound- 
aries and blockages may be located. A single cell 
of the mesh is enlarged in Fig. 3 which shows the 
location of each field variable p, u, v and w rela- 
tive to this (I, J) cell. The pressure and swirl 
velocity y are located at the center of each cell 
and the radial and axial velocities are on the 
boundaries as shown. Thus normal velocities lie 
directly on the physical boundaries of the solu- 
tion domain, while the tangential velocities and 
pressure are displaced half a call interval in- 
*side the flowfield. In this wa*y the exterior fic- 
titious cells are particularly convenient when ap- 
plying the boundary conditions. 

Boundary Conditions 

Finite difference equations simulating the dif- 
ferential equations are set up and solved by way of 
a time-march process applied to cells within the 
flow domain of interest. Cells touching the bound- 
ary thus utilize the value on the boundary (in the 
case of tangential velocities). Interior normal 
velocity calculations take the zero normal wall 
values „ the given normal inlet values, or the yet- 
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Fig. 2 General grid system covering the 
flow domain 





Fig. 3 Arrangement of finite difference 
variables in a typical cell 

to-be-determined outlet values as appropriate during 
their calculations. Interior tangential velocity 
rnlcularinns adjacent to boundaries use the fictiti- 
ous values which are placed in the surrounding layer 
nf complementary cells. Specification of these is 
after each time-step and after each sweep of the 
cells during the pressure-velocity iteration. With 
a coarse grid, fxee-slip conditions are appropriate 
for tangential velocities, and external values are 
set equal to their associated immediately interior 
values. On the other hand, with a fine grid com- 
puting through the boundary layer, no-slip conditions 
are appropriate for tangential velocities and exter- 
nal values are set equal to the negative of their 
associated immediately interior values. In the com- 
putations shown later, free slip boundary conditions 
are applied on all wall boundaries except those cor- 
responding to internal obstacles . These are incor- 
porated via blocking out appropriate cells, setting 
to zero their boundary normal velocities as recommen- 
ded as a simple approach in Kef. 13. If accuracy 
demands, more realistic approaches are easy to incor- 
porate. Specification of all velocities at the out-^ 
let is via the continuative condition. On the axis 
of symmetry the usual zero normal velocity and free- 
slip axial velocity specification are applicable; 
the swirl velocity is given a definite zero value. 


Solution Technique 

Numerical solution of the problem is obtained 
from the time-dependent Navier-Stokes equations, a- 
long with initial and boundary conditions, following 
the Los Alamos SOLA technique. 1 3 Each forward time 


2 



j g%B|g 


if -t .... . ' ’♦ft * * 


step is obtained from the momentim equations, using 
an explicit scheme with donor cell approach for the 
convection terms. The newly-calculated velocities 
will not, in general, satisfy the continuity re- 
quirement* So at each step of the march a pres- 
sure-velocity correction iteration procedure is 
applied; iteration continues with overrelaxation 
until the cell mass imbalance for each cell is 
less than some prescribed small quantity. Criteria 
for the stability and accuracy of the algorithm can 
be found elsewhere , 1 %t l6 “ 18 along with further de- 
tails about the solution technique. Convergence to 
the approximately steady-state solution is estab- 
lished by taking many forward time-steps. 

Sample Computations 
Standard Prediction 

To demonstrate the capabilities of the solution 
scheme, predictions are now made and discussed for 
the case of a blockage located centrally in a sud- 
den expansion chamber, at the expansion plane. In 
general, the blockage ratio BR is (d/D) 2 * 0.25 and 
the chamber expansion ratio D*/D * 1.67. Prechamber 
length A and subsequent main chamber length L are 
0.5 D f and 1,05 D*, respectively. A 12 x 30 square 
mesh system [including fictitious cells] overlays 
the flow domain. A flat inlet axial velocity is 
imposed, with zero radial velocity. Transient com- 
putations are continued until steady state condi- 
tions establish themselves, roughly two flowfield 
resistance times later, corresponding to about 500 
time-steps, and requiring about 120 s of CDC Cyber 
174 CP time. The number of pressure-velocity 
iterations required per time-step depends on initial 
nonuniformities, but typically it is about 200 dur- 
ing early stages of transient computations and de- 
creases to 1 as the steady-state solution is ap- 
proached. 

Predictions are first discussed for the case 
of a blockage located at the expansion plane, with 
sizes as just described. Figures 4 and 5 show the 
predicted recirculation zone size for the case of 
disk and cone blockages. Each computation in these 
and all other results shown is determined at the 
same Reynolds number, so chosen to match only the 
recirculation zone length for flow over the disk 
with blockage ratio BR of 0.25. Thus, the time- 
mean behavior of the turbulent flow is simulated 
via the use of an artificially large constant value 
given to laminar viscosity, in a manner reported 
elsewhere. 17 * 10 Comparison with the experimental 
data of Ref. 10 is used for this single calibration. 
Nevertheless predicted zone lengths for other 
geometric configurations are seen to be very good 
in comparison to the daM. As expected, the 
smaller the blockage forebody geometry cone half- 
angle [in going from disk with a * 90° to a cone 
with a - 45° ] the smaller the recirculation zone 
becomes. Figure 6 shows mass flow recirculated as 
a function of distance downstream for the case of 
a disk blockage with BR * 0.25. Again results com- 
pare favorably with the data. Velocity vectors are 
shown in Fig. 7 where central and corner recircula- 
tion zones are clearly evident. 

Swirl Effect 

The effect of inlet swirl, with flat swirl 
velocity w given by 

W <n ” U <n t3n * 
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Fig* 4 Predicted and measured recirculation zone 
size with disk blockage, nonswirling flow 

[Blockage ratio * 0.25, predictions, 

— experiment 1 0 ] . 





Fig. 5 Predicted and measured recirculation zone 
with 45 deg. cone blockage, nonswirling 
flow [Blockage ratio - 0.25, — — — pre- 
dictions, experiment 10 ]. 



Void 


Fig. 6 Reverse mass-flow rate behind disk, 
nonswirling flow [Blockage ratio * 
0.25. predictions, experi- 

ment 10 ] . 
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Fig. 7 Predicted velocity vectors with disk 
blockage and nonsvirling flow. 

is now considered as a perturbation about the 
standard disk prediction. This assumes that the 
inlet flow passes through swirl vanes so imposing 
a swirl angle <}> to the axis direction. Recircula- 
tion zone size is seen in Fig. 8 to increase as 
swirl vane angle $ increases. This general effect 
of swirl is well-known, although specific experi- 
mental data in the present configuration are not 
available. Corresponding velocity vector plots 
for the <}> ■ 60° case are given in Fig. 9. In 
comparison with Fig. 7 for nonswirling flow, an 
increase in central recirculation zone length and 
a decrease in corner recirculation length are evi- 
dent. These facts are well-known to combustion 
engineers. 1 -*-- In the present case, strong centri- 
fugal effects have also encouraged recirculation 
prior to the disk. Similar swirl predictions, albe- 
it without central blockages, are available using 
more complex turbulent computer codes. 21 * 22 



Fig. 8 Effect of swirl on recirculation zone 

boundary behind disk blockage [Blockage 
ratio - 0.25] 



Fig. 9 Predicted velocity vectors with disk blockage 
and swirling flow with <j> ** 60° 


Blockage Ratio Effect 

Changes in disk diameter are accompanied by 
corresponding inverse changes in inlet axial velo- 
city, so as to maintain the same Reynolds number 
based on disk size. However the velocity passing 
the blockage is then affected by available flow 
area, which is related to the inverse square of disk 
diameter. Prediction of blockage ratio BR * (d/D) 2 
effect on central recirculation zone length [normal- 
ized with respect to particular disk diameter d] is 
given in Fig. 10 for nonswirling /low. Corresponding 
experimental data are very well predicted. 



Fig. 10 Effect of blockage ratio on recirculation 
zone length behind disk blockage 

Lateral Position Effect 


The effect of lateral position of disk blockage 
[with blockage ratio BR - 0.25] is illustrated via 
velocity vector plots in Figs. 11 and 12, which show 
qualitatively the effect of upstream and downstream 
movements, respectively, of the blockage under non- 
swirling conditions. In the former case, higher 
velocities pass the disk and occur near the expansion 
corner, giving rise to larger central and corner 
recirculation zones than those seen in Fig, 7. In 
the latter case, lower velocities with an outward 
radial component pass the disk and occur near the 
expansion corner, giving rise to smaller central and 
corner recirculation zones than those seen in Fig. 

7. Finally, Fig. 13 shows recirculation zone bound- 
aries with respect to the disk location for the three 
axial locations considered. 



Fig. 11 Predicted velocity vectors with upstream 
located disk blockage and nonswirling 
flow 
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Fig. 12 Predicted velocity vectors with downstream 
located disk blockage and nonswirling flow 



Fig. 13 Effect of blockage location on recircula- 
tion zone boundary behind disk blockage 
[Blockage ratio * 0.25] 


Conclusions 

Predictions of recirculation zones behind 
bluff-body flame stabilizers are now possible, 
showing the effects of complex geometric configu- 
rations on the size, shape, and reverse mass flow 
rate in subsequent recirculation zones. These 
recirculation eddies are now amenable to numerical 
solution via the simplified laminar flow finite 
difference technique presented here. Clearly the 
prediction of trends given in this paper take an- 
other step towards providing designers of practical 
equipment with useful qualitative results more 
cheaply and quickly than currently possible by more 
complex analytical and experimental means. 
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Nomenclature 

D Prechamber diameter 

D* Main chamber diameter 

d Blockage diameter 

i, j Mesh point 

Jt P re chanter length 

L Main chamber length 

L o Length of recirculation zone 

p° Ratio of pressure to constant density 

r Radial coordinate 

t Time 

u,v,w Velocity components in x(5r) , y (= axial) 

and 0 coordinate directions 
x,y,0 Cylindrical polar coordinates [in radial, 

axial and swirl directions] 

Y Axial coordinate measured from blockage 

o . 

base 

V Kinematic viscosity 

(f> Swirl vane angle 
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Abstract 


The six-orientation single hot-wire technique 
has been applied to the complex flowfield of a 
swirling, confined jet. This flowfield, which 
contains a rapid expansion with resulting recircu- 
lation regions, is typical of those found In gas 
turbine engines and ramjet combustors. The pre- 
sent study focusses on turbulence measurements in 
such a flowfield In the absence of chemical 
reaction. 

A modification to the six-orientation hot-wire 
technique developed by King has been made, which 
Incorporates the determination of turbulent shear 
stresses (in addition to normal stresses) and 
ensemble averaging of redundant turbulence output 
quantities. With this technique, flowflel d surveys 
have been performed in both swirling and nonswirl- 
ing axi symmetric confined jets. Where Independent 
data exist, comparisons have been made which 
demonstrate the reliability of the technique. 
Finally, a sensitivity analysis of the data reduc- 
tion technique has been completed which forms the 
major ingredient in an uncertainty analysis. 

Nomenclature 


A,B,C Calibration constants In Equation 1 

AO, 130, CO Cooling velocity functions in Tabic 1 
0 Test section diameter 

d Inlet nozzle diameter 

E Hot-wire voltage 

U Velocity function for axial velocity 

W Velocity function for azimuthal velocity 

V Velocity function for radial velocity 

G Pitch factor 

K Yaw factor 

K, , Covariance for cooling velocities 

P 0 Z p , and Z Q 

P,Q,R Selected hot-wire probe positions 

Rej Inlet Reynolds number 

u Axial velocity \ in 

v Radial velocity i coordinates 

w Azimuthal (swirl) velocity )y>n facility 

u.w.v Three components of velocity In probe 

coordinates defined by Figure 5 
x,r,e Axial, radial, azimuthal cylindrical 

polar coordinates 

Z Effective cooling velocity acting on a 

wire 
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a Side-wall expansion angle 

Y z z Correlation coefficient (estimated) 

1 j between cooling velocities of 

adjacent wire orientations 
a 2 Variance of a given quantity 

<fr Inverse function of calibration equation 

$ Swirl vane angle 


Subscripts 

1,2, 3, 4, 5, 6 Refers to the six probe measuring 
positions 

1 J Dummy Indlcles which take the values 

1 to 3 

P,Q,R Refers to the three selected 

measuring positions 

rms Root-mean- squared quantity 

Superscripts 

T1me*mean average 
1 Fluctuating quantity 


1. Introduction 

1.1 The Gas Turbine Combustor Flowfield 


Recent emphasis on fuel economy and pollutant 
suppression has sparked a renewed Interest In gas 
turbine combustor analysis. A typical axi sym- 
metric gas turbine engine combustor Is shown in 
Figure 1. Flowflelds within such combustors typi- 
cally have a rapid expansion and strong swirl Im- 
parted to the Incoming air, which result In corner 
and central recirculation regions. The swirling, 
recirculating, turbulent flows within combustors 
present one of the more difficult fluid dynamic 
problems to analyze. This complexity Is increased 
many fold by the processes of combustion and heat 
transfer within the flowfield. Despite the com- 
plexity of combustor flows, significant progress 
Is being made In their analysis. 1 



Figure 1. Axi symmetric Combustor 
of a Gas Turbine Engine 
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The present paper reports on research which Is 
part of an extensive experimental and computa- 
tional study of gas turbine flowflelds In the 
absence of combustion. Figure 2 shows the charac- 
teristics of the simplified flowfleld being Inves- 
tigated. Flow enters through a jet of diameter d 
Into a tube of diameter D, after being expanded 
through an angle a. Before entering the tube, the 
flow may be swirled by a swlrler located upstream 
of the Inlet plane. Shown schematically are the 
corner recirculation zone (CRZ) and the central 
toroidal recirculation zone (CTRZ) which are typi- 
cally present In these flows. 



Figure 2. Idealized Combustor Flowfleld 


The swirling confined jet flowfleld shown In 
Figure 1 Is being Investigated at Oklahoma State 
University and at Dynamics Technology, with vari- 
ous methods of approach. Analytically, a computer 
program (STARPIC) has been assembled which Is de- 
signed specifically to calculate the swirling con- 
fined jet flowflelds. 2 Experimentally, a series 
of flow visualization experiments coupled with 5- 
hole pitot probe measurements have been used to 
characterize the time-mean flowfleld. 3 Hot-wire 
measurements of the turbulence properties are also 
being conducted. This paper reports on the Ini- 
tial results of the hot-wire measurements In the 
confined jet flowfleld. 

Several studies on time-mean flowflelds of the 
type just described have been carried out using 
various turbulence measuring techniques. ‘ 4_1T) 
Unfortunately, most of the techniques used do not 
give complete and detailed information about the 
flow In terms of all its time-mean and turbulence 
quantities. In addition, no experiments have been 
performed on the specific geometry of the present 
study in the presence of Inlet swirl. To develop 
further the flowfleld computational techniques, 
including the turbulence modeling, there Is a 
strong need to obtain experimental estimates of 
the turbulence and mean flow quantities In such 
flows. 

1.2 The Turbulence Measurement Problem 

Turbulence measurement In a complex flowfleld 
has always been a complicated problem encountered 
by engineers. In the past, turbulence phenomena 
have been discussed by various authors In detail 
and various methods of turbulence measurements 


have been suggested. 11 " 13 One of the most widely 
used Instruments to obtain turbulence quantities 
Is the hot-wire anemometer, the most common of 
which Is a single hot-wire. When used at a single 
orientation and In a two-dimensional flow with a 
predominant flow direction, a single hot-wire can 
measure the streamwlse components of the time-mean 
velocity and the root-mean- square velocity fluc- 
tuation at a particular location In the flowfleld. 
A two-wire probe can be used to determine the 
time-mean velocities, streamwlse and cross stream 
turbulence Intensities, and the cross correlation 
between the two components of the velocity fluc- 
tuatlons . 11 *' 15 

Hot-wire measurements In complex three-dimen- 
sional flowflelds are considerably more difficult 
than In one- or two-dimensional flowflelds In 
which the mean flow Is predominantly In one direc- 
tion. To measure the three velocities and their 
corresponding fluctuating components In a three- 
dimensional flowfleld such as encountered In com- 
bustor simulators, there are two methods that can 
be employed at a point In the flowfleld: 

1) A three-wire probe used with a 
single orientation. 

2) A single- or double-wire probe 
used with multi-orientation. 

The three-wire probe technique permits the 
necessary simultaneous measurements from which 
three Instantaneous velocity components can be 
determined. The appropriate signal processing can 
produce estimates of mean velocity components and 
normal and shear turbulent stresses (such as 

u' 2 and uV ). 

The three-hot-wire probe technique Is signifi- 
cantly more complex than the single wire multi- 
orientation techniques. A multi-dimensional probe 
drive Is required to orient the probe in approxi- 
mately the mean flow direction. Also, sophistica- 
ted signal processing electronics Is required to 
handle the three Instantaneous hot-wire voltages. 
Finally, the three-wire probe typically has less 
spatial resolution In comparison with a single 
wire probe. 

Multi-orientation of a single hot-wire Is a 
novel way to measure the three components of a 
velocity vector and their fluctuating components. 
A method devised by Dvorak arid Syred 16 uses a 
single normal hot-wire oriented at three different 
positions such that the center one is separated by 
45 degrees from the other two. The velocity 
vector at a location Is related to the three 
orthogonal components using pitch and yaw factors 
as defined by Jorgensen. 17 The data are obtained 
in the form of mean and root-mean-square voltages 
at each orientation. However, the measurements 
done with a single wire do not supply all the 
information needed to obtain the turbulence quan- 
tities. Therefore, in addition to a single wire, 
Dvorak and Syred used a cross-wire probe to obtain 
the covariances between the voltages obtained at 
adjacent hot-wire orientations. 

King 18 modified the technique developed by 
Dvorak and Syred. His method calls for a normal 
hot-wire to be oriented through six different po- 
sitions, each orientation separated by 30 degrees 
from the adjacent one. Mean and root-mean-square 
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voltages are measured at each orientation. The 
data reduction Is performed using some assumptions 
regarding the statistical nature of turbulence, 
making It possible to solve for the three time- 
mean velocities, the three normal turbulent 
stresses, and the three turbulent shear stresses. 

1.3 The Scope of the Present Study 

In the present study, the six-orientation 
single normal hot-wire technique Is being employed 
to obtain the turbulence quantities In the combus- 
tor simulator confined jet flowfleld. Measure- 
ments have been carried out for both swirling and 
nonswirling flow with expansion angles of 90 
degrees (sudden expansion) and 45 degrees (gradual 
expansion). Only the 90 degree angle data are 
presented here and the Reynolds number of the 
Inlet flow Is * 5 x 10 1 * which Is comparable with 
aircraft combustor flows (although our experiments 
are performed In nonreacting flows). The data 
reduction procedure extends King's technique to 
obtain turbulent shear and normal stresses using 
six basic response equations representing the six 
orientations of a normal hot-wire positioned In 
the flowfleld. Certain modifications are made In 
the procedure to calculate covariances which are 
an Integral part of the data reduction procedure. 
An uncertainty analysis Is performed on the tech- 
nique which reveals the sensitivity of this tech- 
nique to various Input parameters discussed In the 
later parts of this paper. Some of the turbulence 
quantities obtained are compared with measurements 
performed by Chaturved1 s using a crossed-wlre 
probe In a corresponding flow situation. 


2. Experimental Facility and Instrumentation 
2.1 Idealized Flowfleld 

The facility, designed and built at Oklahoma 
State University, Is a simulation of a typical 
axisymmetrlc combustion chamber of a gas turbine 
engine shown In Figure 1. The schematic of the 
test facility with the idealized flowfleld Is 
shown in Figure 3. Anblent air enters the low- 
speed wind tunnel through a foam air filter. The 
air then flows through an axial flow fan driven by 
a 5 h.p. varldrlve motor. Thus, the flow rate can 
be varied for different test conditions. The flow 
passes through a turbulence management section 
which has two fine-mesh screens, a 12.7 cm length 
of packed straws, and five more fine-mesh screens. 



Figure 3. Schematic of the Experimental 
Facility 


The axl symmetric nozzle was designed to pro- 
duce a minimum adverse pressure gradient on the 
boundary layer to avoid flow unsteadiness asso- 
ciated with local separation regions. The area 
ratio of the cross sections of the turbulence 
management section to that of the nozzle throat Is 
approximately 22.5. The diameter, d, of the 
nozzle throat is approximately 15 cm, 

The test section Is composed of a swirl er 
(optional), an expansion block, and a long plexi- 
glass tube. The expansion block, attached after 
the swlrler, Is a 30 cm diameter annular disk of 
wood. At present, there are three expansion 
blocks, and the appropriate choice gives a * 90, 
70, or 45 degrees. The flow Is expanded Into a 
plexiglass tube of diameter, D, of 30 cm, thus 
giving a diameter expansion ratio (D/d) of 2. The 
test chamber has no film cooling holes or dilution 
air holes, and the chamber wall of the test sec- 
tion Is a constant diameter pipe. 

2.2 Hot-Wire Instrumentation 

Jhe anemometer used for the present study Is 
DISA type 55M01, CTA standard bridge. A normal 
hot-wire probe, DISA type 55P01, Is used In the 
experiments. This probe has two prongs set 
approximately 3 mm apart which support a 5 pm 
diameter wire which Is gold plated near the prongs 
to reduce end effects and strengthen the wire. 
The mean voltage is measured with a Hlckok Digital 
Systems, Model DP100, integrating voltmeter and 
the root-mean-square voltage fluctuation Is mea- 
sured using a Hewlett Packard, Model 400 HR, AC 
vol tmeter. 

The hot-wire Is supported In the facility by a 
traversing mechanism shown schematically In Figure 
4. It consists of a base that Is modified to 
mount on the plexiglass tube of the test section 
at various axial locations. The hot-wire probe Is 
Inserted Into the tube through a rotary vernier 
and the base. The rotary vernier Is attached to a 
slide which can traverse across the flow chamber. 
Thus, It Is possible for the probe to be traversed 
to any radial location at selected downstream 
locations In the flowfleld and to be rotated 
through 180 degrees. 


Travaric Unit 



Figure 4. Hot-Wire Probe Mounted on the 
Test Section 
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2.3 Calibration Nozzle Z 2 . y2 + 52^2 4. *232 (2) 


The hot-wire is calibrated in a small air jet. 
The facility consists of a compressed air line, 
which delivers the desired flow rate through a 
small pressure regulator and a Fisher and Porter 
Model 10A1735A rotameter. The jet housing con- 
sists of an effective flow management section fol- 
lowed by a contoured nozzle with a 3.5 cm diameter 
throat. A rotary table Is used to hold the probe 
While It is being calibrated in three different 
orientations. 

3. Hot-Wire Data Analysis 
3.1 Hot-Wire Response Equations 

The six-orientation hot-wire technique 
requires a single, straight, hot-wire to be cali- 
brated for three different flow directions In 
order to determine the directional sensitivity of 
such a probe. The three directions and three 

typical calibration curves are shown In Figure 5. 
In these relations, tildes signify components of 
the Instantaneous velocity vector In coordinates 
on the probe. Each of the three calibration 

curves is obtained with zero velocity In the other 
two directions. The calibration curves demon- 
strate that the hot-wire Is most efficiently 
cooled when the flow is In the direction of the u 
component, whereas, the wire Is most Inefficiently 
cooled when the flow is in the direction of 
the w component. Each of the calibration curves 
follows a second order, least square fit of the 
form: 

E| * Aj + B i Uj/ Z + CjUj (1) 


which Is an extension of the commonly used King's 
law. In this equation, , Bw, and Cj are 
calibration constants and u, can take on a value 
of u , v, and w for the three 1 calibration curves, 
respectively. 



where G and X are the pitch and yaw factors 
defined by Jorgensen 17 to be: 

G , §021^21 , and 
u( v and vAo) 

K a v(u and w»0) (3) 

w(u and v*0) 

which are evaluated from the three calibration 
curves {Figure 5) for a constant value of E z . 
Equation (3) shows that the pitch and yaw factors 
are calculated with the v component 1 * 2 In equa- 
tion (1) of the effective cooling velocity as the 
reference. Therefore, the calibration constants 
used In equation (1) are the coefficients In the E 
vs. v calibration of Figure 5., I.e., In a general 
flowfleld: 

E 2 = Ag + B 2 Z 1/2 + C 2 Z 

with Z as given In Equation (2) above. 

Figure 6 shows the pitch and yaw factors as a 
function of hot-wire voltage determined from the 
calibration curve of Figure 5. Both factors vary 
with hot-wire voltage, but the yaw factor Is far 
more sensitive. The sensitivity analysis dis- 
cussed In the next section demonstrates that 
uncertainties associated with the varying pitch 
and yaw factors do not seriously affect the accur- 
acy of the estimated flow quantities. 



Figure 6. Pitch and Yaw Factors Plotted Against 
Hot-Wire Mean Effective Voltage 


Figure 5. The Three-Directional Hot-Wire 
Calibration 


When the wire Is placed In a three dimensional 
flowfleld, the effective cooling velocity experi- 
enced by the hot-wire Is: 


To carry out measurements In the confined Jet 
flowfleld, the wire Is aligned In the flow In such 
a way that In the first orientation, the wire Is 
normal to the flow in the axial direction and the 
probe coordinates coincide with the coordinates of 
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the experimental facility. Thus, the six equa- 
tions for the Instantaneous cooling velocities at 
the six orientations, as given by King 18 are: 


probe. Therefore, Equations (S) must be expressed 
In terms of mean and root-mean-square values. 
Equation (1) can be written as: 


Z 2 * v 2 + G 2 u 2 + K 2 w 2 

f 

Z| - v 2 + G 2 (u cos 30°+ w sin 3Q 0 ) 2 ' 

+ K 2 (w cos 30°- u sin 30 0 ) 2 

Z| ■ v 2 + G 2 (u cos 60°+ w sin 60°) 2 
+ K 2 (w cos 60°- u sin 60°) 2 

Z| * v 2 + G 2 w 2 + K 2 u 2 (4) 

Z 2 * v 2 + G 2 (w sin 120°+ u cos 120 0 ) 2 
+ K 2 (u sin 120°- w cos 120°) 2 

Z^ ■ v 2 + G 2 (w sin 150°+ u cos ISO 0 ) 2 
+ K 2 (u sin 150^- w cos 150° ) 2 


z r 4C 2 (v e i^ 1/2/2C 2^ ' (6) 


The above equation is in terms of Instantaneous 
velocity Zf and Instantaneous voltage Ej. In 
order to obtain an expression for time-mean velo- 
city as a function of time-mean voltage, a Taylor 
series expansion of Equation (6) can be carried 
out as follows: 


V z l(V E J) = 4(Ei)+t{-|^-+ 7}-^-+.. 
where $ « Z^(E^). 


(7) 


The Taylor series Is truncated after second order 
terms assuming the higher order terms to be rela- 
tively small. Time averaging both sides of the 
above equation and employing the fact that E'=0, 
yields: 


Solving simultaneously any three adjacent equa- 
tions provide expressions for the Instantaneous 
values of the three velocity components, u, w, and 
v In terms of the equivalent cooling velocities 
(Z,, Z 2 and Z 3 for example, when the first three 
equations are chosen). Thus, the general form of 
the Instantaneous velocity components Is given as: 


♦ + 


1 5 2 <fr 
7 &E 2 



( 8 ) 


To obtain Zl 2 = a? , the relationship as given by 
Hlnze 19 Is: 1 l i 


nn2 1/2 . 1/2 

U - [{A0 + (AO 2 +521 } l —) 

~ (G 2 -K 2 ) 

„ n2 1/2 . 1/2 

W = [{ -A0 + (AO 2 + - 3 -) } ] 

V « [ CO - q2+K - ) - (A3 2 + ^r-) ,Z ] 1/Z 
(G 2 -K 2 ) 


(5) 


The values of AO, BO and CO depend on the set 
of the three equations chosen and are given In 
Table 1 for appropriate equation sets. 


Z i 2 “ a Z = Ex P ec - (Expec [ Zl ]) 2 . (9) 


Using Equation (8) as the basis, Expec [Zp 
and (ExpecCZ^]) 2 can be evaluated and substituted 
Into Equation (8) to get: 


Z 1 2 


a O : 


2 

i~ aEJ a E 1 


-( 


1 & 2 $ 

? 5E 2 


( 10 ) 


Thus, Equations (8) and (10) give the mean and 
variance of effective cooling velocities In terms 
of the mean and variance of the appropriate vol- 
tages. 


TABLE l 


Values of AO, DO, And CO In Various Equation Sets 


Equation 5ct 
P,Q,R Choice 

m 

00 

CO 

U 1 


t-2Z* + 3Z* - Z*} 

(Zj-ZZ,zp 

2, 3, 4 


(-Z2 . 3Z* - 2ZJ) 

(Z z " Z j f Z*) 

3, 4, A 

(Z» - ZZJ , Z*) 


(zj-zj.zp 

4, 5f 6 


l-nl * n] - 1 \) 

(Zj-zz.zp 

S, 5, 1 

(•z* ♦ 1*1 

t-Zj + 3ZJ - 2Z|) 

czz-zz.zp 

6, l, t 

(-ZJ ♦ ZZj - Z») 

C-z* ♦ z|1 

(zz.zj.zj) 


In a 3-dlmenslonal flow, It Is usually desir- 
able to obtain the mean and variance for the Indi- 
vidual velocity components In axial, azimuthal, 
and radial directions, and also their cross cor- 
relations. The procedure to obtain the mean and 
variance of the Individual velocity components Is 
the same as for the effective cooling velocities 
except that u, w and v are functions of three 
random variables and there are extra terms In the 
Taylor expansion to account for the covariances of 
the cooling velocities, Thus, the axial mean 
velocity component ar given by Dvorak and Syred, 16 
and King 18 Is: 


u-u(z p ,z 0 ,z R ) + »^ t | } an 


However, these equations cannot be directly 
used because It Is Impossible to obtain Zj, Zg and 
Z 3 at a single Instant in time with a single wire 


where K z ^ Is the covariance of the cooling velo- 
cities Z^ and Zj and Is defined as: 


5 










Vj 
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( 12 ) 


Identical expressions 
obtained In terms of 
Derivatives 


for w and v can also be 
W and V, respectively. 

are determined 


of the form a^/dZjdZj 
analytically from equations (5) and Table 1. 

Also, the normal stresses are given as: 


i'2 




• a? + 
L \ 


3 3 

} | 5 77 

M*) 0 1 


au au 


z i z j 


tlons would be such that their contribution to the 
cooling of the wire would be related by the cosine 
of the angle between the wires. This assumption 
leads to the following three values of the corre- 
lation coefficients. 


Y 7 7 

Vq 


= cos 30 * 0.9 


(16) 


y 7 7 = cos 30 » 0.9 

Vr 


To relate Yz p z R with Yz p Z Q a™»Vz Q z R . King Intro- 
duced the following relationship: 


r 1 r a 2 l) 

* ill ^ 


3 

+ I 

4 1 1<j 


a 2 u 

■aT^STJ 



(13) 


with similar expressions for w' 2 and v' 2 . 

Finally, the expressions for shear stresses as 
given by Dvorak and Syred 16 are of the form: 


177'= y mm 2 m K 


- [ i 1 *! | • K, z ) 

2 1=1 bZ 2 Z 1 1 1<;j j aZ 1 aZ j Z i Z j 


. [ l | a 2 v . _ 2 A % S a 2 v 


i + f (14) 


1=1 az 2 "1 i^.j 


% Z R “ ” % Z Q Yz Q Z R 

where n Is given a value of 0.8. 
becomes : 

y 7 7 = (0.8)(0.9)(0.9) = 0.65 
4 P 4 R 

The three covariances are then obtained by substi- 
tuting the corresponding values of the correlation 
coefficients into Equation (15). 

The present study, however, uses Equations 
(16) and (18) during the entire data reduction. 
The reason for this is contained In the results of 
the sensitivity analysis presented In the next 
section. This analysis demonstrated that there is 
not significant error magnification In the data 
reduction due to the correlation terms. 


(17) 

Hence vz p z R 

(18) 


4. Results of Hot-Wire Measurements 


Expressions for u ! w* and v'w* can also be obtained 
in a similar manner and are given In Reference 22. 

3.2 Calculation of Covariances 


Dvorak and Syred 16 used a DISA time correlator 
(55A06) to find the correlation coefficients 
between the velocity fluctuations In the three 
directions. King's approach Is to use the Infor- 
mation obtained by all six orientations and devise 
a mathematical procedure to calculate the covari- 
ances. 


Covariances are calculated using the relatlon- 
shl p : 


N z z 

Vj 


Vj 


w. 




1/2 


(15) 


where yz.z* is t ^ 1e correlation coefficient between 
the two cooling velocities and Zj. By defini- 
tion, the absolute value of the correlation coef- 
ficient Y 7 7 is always less than 1. 

King 18 made certain assumptions to calculate 
the covariances. However, he observed that at 
times the calculated value of the correlation 
coefficient Is greater than one at which instance 
he assigned previously fixed values to the corre- 
lation coefficients. He argued that if two wires 
are separated by an angle of 30 degrees, the fluc- 
tuating signals from the wires at the two loca- 


The six-orientation hot-wire technique was 
employed to measure the turbulence quantities for 
swirling and non-swirling flow conditions In the 
confined jet facility described earlier. Also, an 
extensive sensitivity analysis of the data reduc- 
tion was conducted to assist the estimation of the 
uncertainties In the output quantities. 

4.1 Uncertainty Analysis 

The uncertainty analysis Includes a deter- 
mination of the sensitivity of the six-orientation 
hot-wire data reduction to various input param- 
eters which have major contributions In the 
response equations. Pitch and yaw factors (G and 
K) are used In the response equations described in 
Section 3 in order to account for the directional 
sensitivity of the single hot-wire probe. Figure 
6 shows the pitch and the yaw factors plotted 
against the hot-wire mean effective voltage. Both 
the pitch and yaw factors are functions of the 
hot-wire mean effective voltage, but the yaw 
factor is far more sensitive. A one percent 
Increase in the hot-wire voltage reduces the 
pitch factor by 1.3 percent and the yaw factor by 
56 percent. For the present study, the values of 
these factors are chosen at an average hot- 
wire voltage experienced In the flowfield. 
This was a pprop riate since the output quantities 
(u, u^ s , uV, etc) are only weakly dependent on 

the value of K. This can be seen in the data of 
Table 2 which summarizes a sensitivity analysis 
performed on the data reduction program at a 
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representatl ve position in the flowfleld. 


Table 2 demonstrates the percent change In the 
output quantities for a 1 percent change in most 
of the important input quantities. For the data 
presented in this table only quantities ^calculated 

fror. ihe probe orientation combination Zg, Zg and 
are used, for simplicity. In this swirling 
flow Zg was the minimum of the 6 mean effective 

cooling velocities. King 18 has argued that the 
probe orientation combination approximately cen- 
tered around the minimum effective cooling velo- 
city produces more accurate estimates of calcu- 
lated turbulence quantities, than do the other 
orientation combinations. 


tlons of data from adjacent wire orientations. 
One measure of the uncertainty In the output quan- 
tities can be obtained by examining the variance 
In these quantities calculated from the six dif- 
ferent position combinations. Table 3 shows these 
comparison data^for a representative position in 
the flowfleld. For each of the output quanti- 
ties, an ensemble mean x Is calculated together 

with ari_ ensemble standard deviation a. The 

ratio a/x Is a measure of the uncertainty in the 
output quantity. In this table, NR stands for 

not resolved 1 , a problem that occurs when the 

data reduction program attempts to take the square 
root of a negative quantity. In addition, quan- 
tities which are more than three standard devia- 
tions outside the mean are rejected as spurious 
calculations. 


TABLE 2 

Effect cf Input Pirinwtert on Turbulence Qjaniltles 




t CHANGE IN TIMC-HEW AMO TURBULENCE QUANTITIES 

PARAMETER 

t CHANGE IN 
PARAMETER 


7 

w 

W rms 


"rms 

uV 

uV 

v V 


♦1 

♦ 16.10 

*0.66 

*4.98 

♦IS.75 

-2.06 

♦2.75 

♦ 6,0 

♦51.43 

♦11.94 

h 

** 

♦249 

-2.21 

♦11.49 

-6.50 

*2,42 

♦12.80 

♦ 4.0 

♦14,29 

♦7.46 

h 

♦ l 

•10.59 

-0.36 

-0.50 

-1.80 

♦ 7.07 

-9.54 

• 6,0 

-54.29 

-11.94 

£ U 

♦ 1 

* 0 . 2 ; 

-0,06 

♦ 0.14 

♦1.63 

♦0.13 

♦ 0,39 

*2.0 

♦2.86 

♦ 1.49 


♦ 1 

♦ 0.05 

0.0 

♦ 0.14 

0,0 

-0.13 

♦ 1,57 

0.0 

0.0 

♦ 1.49 


♦ I 

-0.16 

*0.10 

-0.14 

-0.63 

. +1.03 

-1,08 

-2,0 

-5.71 

0.0 

5 

♦ l 

-1,0* 

0.0 

-1.01 

-1.0 

0,0 

-0.90 

-2.0 

-2.86 

-1.49 

K 

♦ 1 

♦o.oi 

-0.04 

♦0,01 

*0.01 

0.0 

♦o.oi 

0.0 

0.0 

OrP 

% Z Q 

♦i 

♦ 0.05 

0.0 

♦0.14 

-0.13 

-0.13 

-1,77 

j 0.0 

•2.06 

♦ 1.49 


+i 

*0.21 

1 

♦0.01 

+0.05 

: 

-1.63 

♦0.13 

-049 

0,0 

-*■>71 

!' 

♦ 1.49 

V* 

*1 

j -0.16 

♦0.10 

-0.09 

♦0.13 

o.o 

+0.69 

-2.0 

1 +2.86 

0.0 


It is not unusual in hot-wire anemometry to 
have the mean velocity components and turbulence 
quantities that are measured, be quite sensitive 
to changes in mean hot-wire voltage. For inter- 
pretive purposes, the mean hot-wire voltage varia- 
tions can be thought of as being either errors in 
measuring the mean voltage, or shifts In the indi- 
vidual wire calibrations due to contamination or 
strain ‘aging 1 of the wire. The data of Table 2 
demonstrate that the most serious inaccuracies In 
the measurement and data reduction technique will 
be in the estimates of turbulent shear stresses, 

the most Inaccurate output term being uV. 

♦ 

As already discussed in Section 3, an ad hoc 
assumption is made regarding the numerical values 
of the correlation coefficients used In the deduc- 
tion of time-mean and turbulence quantities. The 
results of the sensitivity analysis (Table 2) show 
the time-mean and turbulence quantities to be 
relatively insensitive to variations in the cor- 
relation coefficients. Therefore, the major ad 
koc assumption made in the technique does not seem 
to have a great effect on the output quantities 
compared to the effect of other input quantities. 

As mentioned earlier, turbulence quantities 
(the output) can be calculated from six combi na- 


TABlE 3 


Setter Anoflj the Turbulence Quantities When Solved by 
Sin Ot f ferent Combinations 


TURBULENCE 

TURBULENCE QUANTITY SOLVED BV SIX COHBlHATtdNS 

MEAN 

Sr/WOMD 


QUANTITY 

1,2.3 

2.3.4 

3,4,5 

4,5,6 

5,6.1 

54+2 

7 

DEVIATION 

a 

o/7 

u/u 0 

0.21 

0.20 

0.21 

0.21 

049 

0.18 

0.20 

0,01 

0.06 

7/7 0 

0 . 1 * 

NR* 

041 

047 

047 

047 

044 

0,04 

0,26 

7/7 

Q.4& 

0,39 

0.39 

0.3S 

0,37 

0,40 

0,33 

0,01 

Q,03 

■W“« 

0,14 

0.14 

0,14 

0.07 

0,08 

0.08 

0,11 

0.03 

0.31 


0.06 

O.H 

041 

0.08 

e,oa 

0,09 

0.09 

0.02 

0.23 


0,13 

0,16 

040 '| 

0.11 

.740 

042 

0,12 

0.02 

0.20 

uV/u* 

NR* 

NR* 

0.012 

MR* 

0.005 

0.004 

0.007 

0,004 

0.62 

—-— 7—2 

u “/% 

0.002 

0.010 

0.002 

0.004 

0.004 

0.000 

0.005 

0.003 

0.72 


0.003 

NR* 

0.003 

0,003 

0.007 

0,001 

0.003 

0,002 

0.58 


* Not Resolved 


The data In Tables 2 and 3 can be used to pro- 
duce estimates In the uncertainties of the calcu- 
lated turbulence quantities. The data suggest 
that uncertainties on the order of 5 percent are 
to be expected In the mean velocity component 
estimates. Normal turbulent stress estimates 
( u 1 rms , etc.) have uncertainties on the order of 
20 to 30 percent and turbulent shear stress esti- 
mates are significantly higher, although most of 
this Is a consequence of taking a product of terms 
such as u‘ and v’ . 

These uncertainty estimates are considered to 
be somewhat conservative. More accurate estimates 
are quite difficult to obtain because, to our 
knowledge,, similar measurements have not been per- 
formed with any other Instrumentation system in 
this geometry flowfleld. Also, comparisons of 
several representative points with Independent 
measurements suggest that the ensemble averages 
estimates are typically in closer agreement than 
are selected sets of three orientations. There- 
fore, all turbulence estimates presented in this 
paper are calculated from ensemble averages of six 
groups of three adjacent wire orientations. Any 
data not resolved are not Included In this avera- 
ging. This approach represents a departure from 
the technique developed by King 18 who typically 


* x/D-1.5, r/D=0.25, *=38 deg (swirling flow). 
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selected one group of three orientations from 
which to calculate his turbulence estimates. 


4.2 Results of Flowfleld Surveys 

Radial distributions of time-mean velocities, 
turbulent normal stresses and shear stresses are 
obtained for both nonswirling and swirling con- 
ditions, at various axial locations In the flow- 
field. 

Nonswlrllng Fl ow. In the confined jet, the 
experiments have been conducted with expansion 
angles of 90 degrees (sudden expansion) and 45 
degrees (gradual expansion) and the results for 
both cases are presented In Reference 22. In the 
Interest of brevity, only the data for a 90 degree 
expansion are presented here. 

Figure 7 shows the radial distribution of 
time-mean axial and radial velocity components at 
various axial locations. The axial velocity dis- 
tributions are compared with a similar study per- 
formed by Chaturvedl 5 with a crossed hot-wire 
probe. Because of the Inability of the six-orien- 
tation hot-wire technique to determine the sense 
of the flow direction In a nonswlrllng flow, the 
presence of the corner recirculation zone was 
observed hy a sudden Increase In the axial velo- 
city closer to the wall. 


Figur.. 0 shows the radial distribution of 
axial and radial components of the turbulence 
intensity at various axial locations In the con- 
fined jet flowfleld. These turbulence Intensity 
components are compared with Chaturvedl 's measure- 
ments 5 and reasonable agreement is found. In 
fact, the agreement In most cases Is better than 
the uncertainty estimates derived from the data 
reduction sensitivity analysis. 

Included on Figure 8 are measured turbulent 
shear stress component (u'v'/U 3 ) profiles for the 

nonswirling confined jet. For the most part, 
these measurements are In reasonable agreement 
with those made by Chaturvedl 5 with a crossed wire 
probe. The two significant exceptions to the good 
agreement occur at the furthest upstream and fur- 
thest downstream locations. Upstream, at x/D - 
0.5 the shear layer Is very thin and, therefore, 
matching data from several wire orientations 
obtained at somewhat different times may be prac- 
tically difficult. We believe the overly large 
measured turbulent shear stress on the centerline 
at the furthest downstream station (x/D = 3.0) to 
be a consequence of the transient nature of the 
flow. The recirculation regions In the confined 
jet oscillate somewhat at a low frequency, likely 
characteristic of the main acoustic modes In the 
tube. These large scale oscillations can have 
significant correlated velocity fluctuations (such 

as uV). 

Swirling Flow The measurements performed In 
the swirling flow are with a = 90°, $ = 38°, and 
x/D * 0.5, 1.0, and 1.5. The object of these 
limited number of experiments was to evaluate the 
reliability and accuracy of the six-orientation 
hot-wire technique before making extensive use of 
the technique. 


The hot-wire results In the case of time-mean 
axial and azimuthal (swirl) velocities shown In 
Figures 9 and 10, are compared with five-hole 
pitot probe measurements performed by Rhode In 
the same experimental facility. Agreement among 
ohe two studies Is fairly good. King 18 suggested 
a method to determine the sense of the axial velo- 
city. He advised comparing the magnitudes of 23 
and Zc given by Equations 4. In the present flow- 
field; the swirl velocity Is always positive and 
the two equations giving Z 3 and Zg differ only In 
the sign of axial velocity. Therefore, when Z 5 
is great r than Z 3 , the axial velocity Is nega- 
tive, otherwise It Is positive. With the proper 
sense being assigned to the x velocity mean com- 
ponent, the presence of central torlodal recir- 
culation zone Is evident in the results of both 
measurement techniques. 

Figure 11 shows the radial distribution of the • 
time-mean radial velocity at various axial loca- 
tions for a swirl vane angle of 38 degrees and 
iii.11 expansion angle of 90 degrees. Flow visuali- 
zation and five-hole pitot probe measurements per- 
formed In Rhode's study 3 show the time-mean radial 
velocity to be negative at axial locations greater 
than x/D = 0.5. In spite of the Inability of the 
six-orientation hot-wire technique to determine 
the sense of the radial velocity, the data are 
presented with the appropriate sign change. There 

Is a reasonable agreement among the two studies In 
measurements of time mean radial velocities except 
at the Initial measurement station. 

Figure 12 shows the radial distribution of 
axial, radial and azimuthal turbulent intensities 
at three axial locations presented. At axial 
locations closer to the Inlet of the confined jet, 
the axial turbulence Intensity is fairly high, up 
to 32 percent for x/D = 0.5 which Is due to the 
large axial velocity gradients closer to the wall. 
However, In the case of radial turbulence inten- 
sity, the profiles are rather flat. The mean azi- 
muthal velocity also experiences sudden changes in 
gradients and, hence, the outcome Is a large azi- 
muthal turbulence Intensity closer to the wall at 
x/D = 0.5. 


Figure 13 shows the shear stresses u'v'/U. , 

2 ... .... . - 2 0 

u'w' /U , and v'w' /U as a function of radial and 

axial distance. The sensitivity analysis showed 
that we should expect large uncertainties asso- 
ciated with evaluation of turbulent shear stresses 
using the six-orientation technique. Therefore, 
the reliability of the profiles of these shear 
stresses shown in Figure 13 Is uncertain at 

this time. Nevertheless, stresses u'v'/U and 

2 0 

u'w'/U 0 are found to have large values closer to 

the wall that one would expect due to steep axial 
and azimuthal velocity gradients. The fact that 
we have found no other measurements of this type 
In a swirling, recirculating flow attests to the 
fact that accurate measurements In such a flow are 
quite difficult. 

Cl osure 

The six-orientation hot-wire technique Is a 
relatively new method to measure time-mean velo- 
city components and turbulence quantities In com- 
plex three-dimensional flowflelds. Applied In 
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Radial Distributions of Time-Mean Axial 
and Radial Velocity Components In the 
Non-Swirling Confined Jet; (Note the 
Difference in Velocity Scales) 
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(5-hole pitot probe) 
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Figure 9. Radial Distribution of Time-Mean Axial 
Velocity In Swirling Confined Jet for 
a = 90° and « * 38° 
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Figure 10. Radial Distribution of Time-Mean 
Azimuthal Velocity in Swirling 
Confined Jet for o * 90° and 4> ■ 38® 


this study to nonreacting axl symmetric flowflelds, 
measurements of time-mean and root-mean-square 
voltages at six different orientations contain 
enough information to obtain the time-mean velo*- 
cities, turbulence intensities and shear stresses. 
At each location In the flow, there are six dif- 
ferent values of each of the above quantities that 
can be obtained using six sets of measurements of 
three adjacent orientations. Ensemble averages of 
the output quantities from the six combinations of 
data appear to produce estimates with the best 
agreement with independent measurements. 


FI owfield surveys of both swirling and non- 
swirling confined jets have been made with the 
six-orientation single hot-wire technique. These 
measurements have been used to calculate estimates 
of the mean velocity components and the normal and 
shear turbulent stresses. Where independent data 
exist, comparisons have been made which demon- 
strate the reliability of the technique. 

In addition, a sensitivity analysis of the 
data reduction technique has been conducted which 
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forms the major Ingredient In the uncertainty 
analysis* It Is demonstrated that the largest 
uncertainties are to be expected In the turbulent 
shear stress estimates. Nevertheless, In non- 
swirling flows the measured shear stresses are In 
reasonable agreement with previous measurements 
made with a crossed-wlre probe. In swirling flow, 
previous similar measurements have not been found. 
Consequently, the universal accuracy of the mea- 


sured technique applied to swirl flows 1$ still an 
open question. 
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Figure 11. Radial Distribution of Time-Mean 

Radial Velocity In Swirling Confined 
Jet for a « 90* and <X> * 38° 


4 <0 
4* 

iea 


1 V.« tn' 

U rms^ u n 


m • 
« • 


13 *10 ■ 04 
© *fO ■ 1.0 
m/0 • 14 




4 . 1 * 




m « 

O** 

014 


tTt - 


1 Alt 4.2*' Lit 





Am 


An 2 TT* Au 

v*w'/u! * 10" 1 

o 


1.0 


Figure 12. Radial Distribution of Turbulent Shear Stresses uV/u\ uV /u* 

2 ® 0 * 

and vV/u Q In Swirling Confined Jet for a = 90° and ® = 38° 



oRici.m pass y 

OF POOR QUALITY 


References 

LI Hey, D.G., "Flowfleld Modeling In Practical 
Combustors: A Review," Journal of Energy, 3, 
No. 4, 1979. “• “ 

L1 1 ley, D.G. and Rhode, D.L., "A Computer Code 
for Swirling Turbulent Axlsymnetrlc Recircula- 
ting Flows In Practical Isothermal Combustor 
Geometries," NASA Contractor Report 3442, 
February 1982. 

Rhode, O.L., LI lley, D.G. and McLaughlin, 
D.K., "Mean Flowflelds In Axlsynmetrlc Com- 
bustor Geometries with Swirl," AIAA Paper No. 
82-0177, 1982, AIAA Journal (in press). 

Krall, K.M. and Sparrow, E.M., "Turbulent Heat 
Transfer In the Separated, Reattached, and 
Redevelopment Regions of a Circular Tube," 
Journ al of Heat Trans fer, pp. 131-136, 
February 1966. 

Chaturvedi , M.C., "Flow Characteristics of 

Axlsymnetrlc Expansions," Proceedings, Journal 
of the Hydraulic Division, ASCE, 69, No. HY3, 
pp. 61-92, 1963. 

Phaneuf, J.T, and Netzer, D.W., Flow Charac- 
teHstlcsIn Solid Fuel Ramjets^ Report No . 
NPS-B7Nt740Sl. Prepared for the Naval Weapons 
Center by the Naval Postgraduate School, 
Monterey, California, July 1974. 

Back, L.H. and Roschke, E.J., "Shear Layer 
Flow Regimes and Wave Instabilities and 
Reattachment Lengths Downstream of an Abrupt 
Circular Channel Expansion," Journal of 
Applied Mechanics pp. 677-681, September 1972. 

Roschke, E.J. and Back, L.H., "The Influence 
of Upstream Conditions on Flow Reattachment 
Lengths Downstream of an Abrupt Circular 
Channel Expansion," Journal of Biomechanics, 
9, PP. 481-483, 1976. 

Ha Mlnh, H. and Chassalng, P., "Perturbations 
of Turbulent Pipe Flow," Proceedings, Sympo- 
sium on Turbulent Shear Flows , Pennsylvania 
State University, pp. 13. 9-13.17, April 1977. 

Moon, L.F. and Rudinger, G., "Velocity Distri- 
bution In an Abruptly Expanding Circular 
Duct," Journal of Fluids Engineering, pp. 226- 
230, March'1977: 


11 Bradshaw, P., An Introduction to Turbulence 
and Its Measurement, Pegamon Fress, New York, 
T97TI 

12 Beer, J.M. and Chlgler, N.A. Combustion Aero- 
dynamlcs, Halsted Press Division, John Wiley S 
Sony Inc., New York, 1972. 

13 Syred, N., Beer, J.M. and Chlgler, N.A., "Tur- 

bulence Measurements In Swirling Recirculating 
Flows," Proceedings, Salford Symposium on 
Internal Flows , Condon, England: Inst. oT 

Mechanical Engineering, pp. B27-B36, 1971. 

llf Wygnanskl, I, and Fiedler, 11., "Some Measure- 
ments In the Self Preserving Jet." Journal of 
Fluid Mechanics , 38, p. 577, 1969. 

15 Pratte, B.D. and Keffer, J.R., "The Swirling 
Turbulent Jet," Journal of Basic Engineering, 
94, pp. 739-748, December' 1972. 

16 Dvorak, K. and Syred, N., "The Statistical 
Analysis of Hot Wire Anemometer Signals In 
Complex Flow Fields," DISA Conference, Univer- 
sity of Leicester, 1972. 

17 Jorgensen, F.E., "Directional Sensitivity of 
Wire and Fiber Film Probes," DISA Information 
No. 11, Franklin Lakes, N.J., pp. 31-37, May 
1971. 

18 King, C.F., "Some Studies of Vortex Devices - 
Vortex Anpllfler Performance Behavior," Ph.D. 
Thesis, University College of Wales, Cardiff, 
Wales, 1978. 

19 Hlnze, J.O. Turbulence , 2nd Edition, McGraw- 
Hill, New York, 1975. 

20 Habib, M. A. and Whltelaw, J.H., "Velocity 
Characteristics of Confined Coaxial Jets With 
and Without Swirl," ASME Paper 797-WA/FE-21, 
New York, N.Y., December 2-7, 1979. 

21 Gupta, A.K. and Lllley, D.G. , FI owfleld Model- 
Ing and Diagnostics, Abacus Press, Yunbrldqe 
Well s', England, 1992 (In press). 

22 Janjua, S.I., "Turbulence Measurements In a 
Complex Flowfleld Using a Six-Orientation Hot- 
Wire Probe Technique," M.S. Thesis, Oklahoma 
State University, Stillwater, Oklahoma, 
December 1981. 


